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ABSTRACT 


Electrons spiraling along magnetic field lines emit electromagnetic 
radiation at radio wavelengths via the cyclotron and synchrotron emis- 
sion processes. Consequently, with the discovery of the Van Allen 
radiation belts, the possibility of detecting synchrotron radiation 
from superthermal electrons trapped in these regions has become an 
important topic of research. Studies, which to the present have been 
limited to frequencies above 2 MHz, show that the radio noise from these 
electrons is not easily detected as it is less than the noise generated 
by galactic and extragalact ic sources. 

In this study, the problem of detecting cyclotron and synchrotron 
noise from these superthermal electrons is analyzed for the frequency 
range 30 kHz to 300 kHz. Due to the earth's ionosphere, ground based 
observation of this noise is improbable. Therefore, the calculations 
are made for an observer in the interplanetary medium. In particular, 
the location is chosen in the geomagnetic equatorial plane at a geo- 
centric distance of 32 earth radii. This position of the observer 
allows the theoretical results to be compared directly with data 
obtained from the Goddard Space Flight Center (GSFC) radio astronomy 
experiment aboard the IMP-6 spacecraft. 

The energy, radiated by the trapped electrons and propagated to 32 
earth radii, is calculated using the equation of radiative transfer. 

To solve this equation, models for the thermal and superthermal electron 
densities, throughout the earth's magnetosphere, are derived from in 
situ and indirect measurements made of these particle populations. 
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The analysis shows that appreciable radio noise in the LF range is 
emitted by electrons in the earth’s outer radiation belt. The inner 
belt does not contribute power at these frequencies since it is located 
in the plasmasphere , a region of high thermal electron density, which 
inhibits radiation at these low frequencies. The intensity of the radio 
noise from the outer radiation belt is positively correlated with K 

P 

(the planetary three-hour- range index). This correlation is due to the 

enhancement of the stably trapped electrons during periods of sustained 

geomagnetic activity. For a frequency f in Hz (30 kHz <i f s: 150 kHz), 

the average calculated brightness, 1(f), is given by a constant, M, 

times the frequency raised to the (-a) power, where both M and (a) vary 

with K as follows: 1 <, K £2, log (M) = -12.82, a = 1.66; 

P P 

3 ^ K ^ 4 , log (M) = -11.50, a = 1.71; 5 £ K £ 6 + , log (M) = -10.13 

" p 

— 2—1 —1 

a = 1.82. The units of 1(f) are watts m Hz steradian 

These theoretical results were compared with data from the GSFC 
experiment aboard the IMP-6 , which detected radio noise in the LF range 
emanating from the earth. The predicted increases in M and (a) with K 

P 

were found in the experimental data. Furthermore, the calculated 

increases in these quantities are within 10% of those seen in the IMP^ 

data, and the experimental and theoretical values of (a) are within 9% 

for each K range. 

P 

The theoretical model, however, underestimates the noise intensity 
observed on the IMP-6 by a consistent factor of 5. This factor is felt 
to be within modeling and experimental errors. Consequently, it is 
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concluded that detectable cyclotron- synchrotron noise in the LF range is 
generated by electrons trapped in the earth* s outer radiation belt, and 
the intensity of this noise increases with increasing geomagnetic 
activity . 

Based on this conclusion, three applications of the results of this 
study are suggested. First, the radio noise can be used to determine 
the spatially integrated flux spectrum for the electrons trapped in the 
earth T s outer radiation belt. Second, the noise can be used to obtain 
an estimate of the maximum interplanetary electron density between the 
observer and the earth. (An example of this application is given in the 
study.) Finally, the ray path calculations in the analysis show that 
a radio occultation experiment between two satellites can be used to 
determine the location of the plasmapause as a function of and/or 
time . 
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UNITS 


Aside from the exceptions noted here, rationalized MKS units are 
employed throughout the study. In the field of geomagnetism it is com- 
mon to refer to the geomagnetic field intensity (H) and flux density 
(B) in gaussian units of oersteds and gauss, respectively. Consequently, 
B is expressed in gauss to conform to common usage in the literature. 

Similarly, in electron physics, electron energies are commonly given in 

3 

units of electron volts (ev), kilo electron volts (kev = 10 ev) and 

g 

mega electron volts (Mev = 10 ev) . Thus electron energies are referred 
to in electron volts rather than in joules. It will also prove conven- 
ient to measure large distances in units of earth radii. To make con- 

-4 2 -19 

versions: 1 gauss = 10 weber/meter ; 1 electron volt = 1.60 X 10 

0 

joules; and 1 earth radius = 6-38 X 10 meters. In all formulae MKS 
units are used unless specifically noted to the contrary. 
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CHAPTER I 


GYRO-SYNCHROTRON RADIATION MODELS 

A . Introduct ion 

Synchrotron emission from relativistic electrons has been postulated 
as the mechanism responsible for radio noise detected from stars 
[Zheleznyakov , 1967; Thorne, 1962], from the planets of our own system 
[Fung, 1966; Chang and Davis, 1962], and from ext ragalactic sources 
[Zheleznyakov, 1967]. The works of Dyce and Nakada [1959] and Vesecky 
[1967] explored the possibility of detecting synchrotron radiation at 
frequencies above 1 MHz from electrons trapped in the earth 1 s inner 
radiation belt. Dyce and Nakada concluded that, for frequencies 
around 30 MHz, the radiation from these electrons could not be detected 
at the earth due to interference from cosmic sources whose emissions 
are one hundred times as strong. Vesecky found that under suitable 
conditions it might be possible for an observer situated at a geocentric 
distance of 4 earth radii in the geomagnetic equatorial plane, to 
detect synchrotron radio noise between 1 MHz and 10 MHz. His study 
showed that, for frequencies below 1 MHz, the thermal electron popula- 
tion decreased the synchrotron radiation below detection levels, 
while the power law energy spectrum of the trapped electrons in the 
inner belt led to a rapid decrease in the energy radiated at frequencies 
above 10 MHz. 

These studies, however, do not consider the possibility of detecting 
radio noise generated at frequencies below 1 MHz by superthermal 
electrons trapped in the outer radiation belt. As will be shown in 
Chapter II, this portion of the trapped energetic electron population 
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is located in a region, referred to as the plasmatrough , which is 
characterized by its very tenuous thermal electron population. Due 
to this low density, ambient plasma, appreciable power can be 
generated and propagated in the LF range by a combined cyclotron- 
synchrotron mechanism (henceforth referred to as the gyro- synchrotron 
process) . 

Even though radio noise at frequencies below 1 MHz cannot 
penetrate the ionosphere, and hence is not detectable by ground based 
stations, there are an increasing number of radio astronomy experiments 
aboard earth orbiting satellites with the capability of detecting 
radio emissions down to frequencies of tens of kHz. The analysis 
to follow is therefore intended to extend the results of previous 
studies by determining if the radiation from energetic electrons in 
the outer radiation belt can be detected at frequencies below 300 
kHz by these spacecraft. The upper bound of 300 kHz is arbitrarily 
set as a means of differentiating between the investigations pursued 
by previous authors and the one presented here. 

To put the problem into more definite terms, we shall determine 
the intensity and spectrum of the noise between 30 kHz and 300 kHz 
(during both geomagnet ically quiet and disturbed times) that could be 
received by a spacecraft located in the geomagnetic equatorial plane 
at a geocentric distance of 32 earth radii. The choice of observer 
location and the lower bound on the frequency range is not arbitrary; 
it is intended to allow a direct comparison between the theoretical 
results of this paper and the data collected by the Goddard Space 
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Flight Center (GSFC) radio astronomy experiments aboard the IMP-6 
spacecraft, (Data from this satellite, kindly supplied through 
Dr. R. G. Stone and Dr. L. Brown of GFSC, are presented and discussed 
in Chapter IV.) The analysis carried out in this paper is general, 
however, and can be extended to data obtained by other spacecraft 
whose orbits and/or radio astronomy experiments are not the same as 
those of the IMP-6 . 

B. Initial Assumpt ions 

In this chapter we shall examine in detail the process by which 
radio frequency electromagnetic waves are radiated by an ensemble of 
superthermal electrons (immersed in a magnetoactive plasma) via the 
gyro-synchrotron mechanism. To model this problem in these general 
terms, however, would prove to be mathematically intractable due to 
the complicated phenomena related to the generation and propagation 
of electromagnetic waves in a magnetized plasma composed of thermally 
excited electrons, ions and neutral particles [see, for example, 
Truelsen, 1971]. Fortunately, the problem can be simplified on the 
basis of the lowest frequency monitored by the experiments aboard the 
IMP-6 spacecraft and the location in the magnetosphere from where the 
observed radio noise is shown to emanate. These two conditions led 
to the results listed below, which shall be referred to by number in 
the subsequent text . 
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1. We shall prove that the regions in which appreciable noise 
is generated are located in the magnetosphere at distances greater 
than 2,1 earth radii (13,440 km) above the earth’s surface. At these 
heights, the magnetoplasma is composed primarily of hydrogen ions (H + ) 
and electrons [Chapman, 1959], which are assumed to be present in 
equal concentration to preserve space charge neutrality. 

2. As shown in Chapter III, the region of emission not only 

lies at a distance of at least 13,000 km from the earth, but is also 
located between -40° £ £ 40°. (Here \ represents geomagnetic 

latitude referenced to a coordinate system based on a centered dipole 
approximation to the geomagnetic field.) This segment of the magneto- 
sphere is that which is swept out by the shaded portion of Figure 1.1 
when it is rotated 360° about its vertical axis. 

The validity of the centered dipole approximation and the assump- 
tions used in obtaining the bounds for the emission region are 
subjects considered in Chapter II when the magnetospher ic model is 
developed. At this point it should be noted that the maximum electron 
and ion (nonrelat ivistic) gyrof requencies (f 0 and f ± ) encountered by 
an electromagnetic wave propagating in this region are 25 kHz and 
13,9 Hz respectively. Since the lowest frequency received by the 
IMP-6 (30 kHz) is many orders of magnitude above f ^ , the ions can 
be assumed stationary and their effect on the generation and propaga- 
tion of the radio noise can be ignored [Stix, 1962]. The thermal 
plasma can therefore be considered as composed of only one constituent: 
electrons . 
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3. For 30 kHz to 300 kHz electromagnetic waves propagating in 
the segment of the magnetosphere described above, the frequency of 
the radio noise will be above the local upper hybrid frequency (f^) 
of the medium. Resonant conditions, where either of the complex 
refractive indices (|j, + ) tend to infinity, therefore do not occur, and 
the real part of the indices for the ordinary (M* + ) and extraordinary 

waves [Ratcliffe's nomenclature, 1962] will be between zero and 
unity, as shown by Figure 1.2. In this and the subsequent text, the 
real parts of the complex refractive indices will be referred to 
simply as the refractive indices, and they are given by 

n = R (u ) (1.1) 

+ ^ + 

n = — R (K ) (1.2) 

+ «) E + 

where K is the propagation vector for the characteristic waves of 
+ 

the medium and R^(x) denotes the real part of x. 

Furthermore, the energy density of the electromagnetic waves 

-12 

generated in the magnetosphere is small, being less than 10 watts 
per square meter, resulting in a magnetic flux density associated 
with the propagating wave of much smaller magnitude than that which 
exists in the emission region. These conditions permit the medium 
to be approximated as linear [Tanenbaum, 1967]. Consequently, the 
linearized theory of propagation [Ginzburg, 1961] can be used. 

4. St ix [1962] shows that a cold- plasma approximation, in 
conjunction with the linearized theory of propagation, is justified 
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Figure 1.2 COMPLEX REFRACTIVE INDICES FOR A MEDIUM IN WHICH Y< 1 
[From Ratcliffe, 1962]. 
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provided the following equations hold. 


"CD + JW e 

2 

m 

e 

. VVL 


2X T 
L e -i 


» 1 j - . ..,- 1 , 0 , 1 ,. .. 

( 1 . 3 ) 



2 


XT «1 
e 


(1.4) 


where 

0 ) - radian frequency of the electromagnetic radiation 
X = Boltzmann's constant 

= temperature of the electrons in degrees Kelvin 
cw = radian gyrof requency of the electrons 

|R e (K 2 )| = magnitude of the real part of the wave propagation 
vector in the z direction 

|Rg(Ry)| = magnitude of the real part of the wave propagation 

vector in the y direction 

m = rest mass of an electron 
e 


Using Equation (1.2) } an average electron temperature in the 
magnetosphere of = 1300 K [Bauer, 1965], and the geometry shown 
in Figures 1.3 and 1.4, the above inequalities can be recast as 





2 2 

n + cos 0 2XT g 


m e c 2 


3 = .... - 1 , 0 , 1 , . . . 


(1.5) 
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and 



( 1 . 6 ) 


Taking the worst case conditions (n^ =1.0 and 0 = 90°) for the 
quantities on the right hand side of these inequalities gives 




» 6 .6 x 10 


j = 


■ 1 . 0,1 


(1.7) 


1.51 x lo 3 


( 1 . 8 ) 



2 2 

Y = Cl - 0 ) 

P 2 = (v* + v 3 )/c 2 

h = V c 

0 2 = v 2 /c 

U ) = uu /v 

g e T 

r ! = V“g 

r (t) = yr cos (<u t) 
x 1 g 

r (t) = Y r , sin (<u t) 
y 1 g 

r (t) = v (t) 

Z 2 


v (t) = dr (t)/dt 
e e 


Figure 1.3 GEOMETRY FOR THE POSITION AND VELOCITY OF A CHARGED 

PARTICLE SPIRALING ALONG A STATIC MAGNETIC FIELD LINE. 
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Figure 1.4 GEOMETRY FOR THE RADIATION EMITTED BY THE CHARGED PARTICLE 
OF FIGURE 1.3. 


In the portion of the magnetosphere which affects the generation 

and propagation of the radio noise, we find 1 kHz ^ f ^25 kHz. 

e 

Thus 1.2 ^ ^ 300 and the inequalities are always satisfied. 

Therefore, the cold plasma approximation is valid for this study. 

5. For electromagnetic waves of frequency (f) propagating in a 
plasma, where (f) is greater than the upper hybrid frequency of the 
thermal plasma, the medium may also be treated as collisionless 
provided [Ginzburg, I960] 


2 2 
U) « (cu - 


(1.9) 
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where 


00 = radian collision frequency between electrons and hydrogen ions. 

The maximum value for u) encountered in this study is 10 radians 

c 

3 

per second [Liehmon, 1965], while (jd - > 3.14 x 10 radians per 

second. Therefore, the medium can be taken as collisionless, and 
the absorption of electromagnetic waves propagating in this medium 
(due to collisions between electrons and ions) can be ignored. 
Equivalently, the dielectric tensor (e) describing the medium is 
Hermetian [Yeh and Liu, 1972] and the propagation vector is 

either real or purely imaginary (for real, positive u)) [Bekefi, 1966], 
Thus, the wave propagates with no loss of energy due to collisions, 
or it is evanescent and does not propagate at all. 

6. The magnetic permeability of the magnetosphere is taken to 
be that of free-space (^ q ) , which is a good approximation for radio 
waves propagating in this region of the earth T s environment [Yeh and 
Liu, 1972]. 

7 . Gyro- synchrotron radiation causes an electron gyrating in a 
magnetic field to lose energy at an average fractional rate per revolu- 
tion of [Bekefi, 1966] 



P 2 - 0 (1.10) 
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where 


c = speed of light in free-space 
e Q = permittivity of free space 

AW = energy lost by the particle in the form of radiation 
Eyt = total energy of the particle (rest mass energy + kinetic energy) 
= velocity component of the particle along the field line as 
defined in Figure 1.4 

The maximum energy, , considered in this study is of the order 

4 -1 

of 5 Mev, and as stated previously, u) £ 15.7 x 10 rad sec 

— 1 6 

Substituting into Equation (1.10) yields AW/E^, « 6 x 10 . Further- 
more, for a fixed this value decreases as p increases (as will 

be substantiated later in this chapter) . This result implies that 
the particle energy (as well as p and p ) is constant over many 
revolutions; consequently, the gyro- synchrotron model can be used 
with little error. 

8. The lowest frequency received by the IMP-6 radio experiments 
is 30 kHz, which gives a free-space wavelength for these signals of 
10 km. As shown in Chapter II, the region in which the electro- 
magnetic waves propagate is characterized by refractive indices 

bounded by 0.1 ^ n £ 1. Therefore, the wavelength in the medium 
■h ' 

for the 30 kHz signal is less than or equal to 100 km. 

In the models developed in Chapter II for the medium, the magneto- 
sphere is represented by concentric shells whose distances from the 
earth in the geomagnetic equatorial plane differ by increments of 
1000 km. That is, spatial variations in magnetospheric parameters 
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of importance to this study can be taken as constant over distances 
less than this value without affecting the final results. Thus, as 
the scale lengths of variations in these parameters are an order of 
magnitude or larger than the wavelength of the radiation propagating 
in these regions, the magnetosphere can be considered as a slowly 
varying medium. Therefore, the principles of geometric optics can 
be used to investigate the transfer of radiation from the source 
regions to the observer [Bekefi, 1966; Ginzburg, I960]. 

9. Finally, the time variations associated with the properties 
of the magnetosphere will generally be much longer than the periods 
of the signals investigated. As a result, the plasma and the emission 
processes can be assumed time stationary. This does not imply, however, 
that gross changes in the magnetosphere due to geomagnetic activity 
will be excluded from this study, as these changes occur over periods 
of hours to days . 

The above results were derived solely on the basis of the frequency 
range considered in this study and the location in the magnetosphere 
from which the major portion of the radiation will be shown to 
emanate. The importance of this fact is that a much more tractable 
model for the generation of radio noise can be developed in the 
next sections. 
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C . An Overall View of the Problem 


Paragraphs eight and nine of the preceding section reduce the 
problem of calculating the signal generated in the magnetosphere and 
propagated to the spacecraft to one of radiative transfer in geometric 
optics. By this analogy, the electromagnetic energy can be envisioned 
as transported along bundles of curves or rays. For points along a 
ray of this type (Figure 1.5), a quantity known as brightness (I) 
satisfies the following differential equation [Bekefi, 1966] 


d 

d A 


lO) 

n r a) 


+ 


JO) 

n 2 r «) 


K(j» 

n 2 a ) 
r 


KO 


( 1 . 11 ) 


where 

I = brightness of the source at a given point on the ray path 

-2 -1 -1 
(watts m Hz ster ) 

K(jg, ) = absorption coefficient as a function of path length 
(nepers m *) 

J(£) = emission coefficient as a function of path length 
(watts m Hz ster ) 

djg = element of path length along the ray 

n r = ray refractive index which defines the angular variation of 
the energy flux in an anisotropic medium and is given as a function 
of n + by Equation (1.12) below [Bekefi, 1966]. 
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by defining a quantity called the optical depth (t) as 

h 

-ra 1 ) = J ku) dt 


( 1 . 12 ) 


(1.13) 


Equation (1.11) can be solved to yield 


I(o) 

% (o) 

+ 


I (°o) e 


■T («) 


n 2 («) 
r + 


f £212 ill 
j n U) 


J Ct) dl 


(1.14) 


For an observer located in the interplanetary medium and monitoring 
the frequency range 30 kHz £ f ^ 300 kHz, we find 


n (o) « 1 
r + 


n" («>) » 1 
r + 


I(co) = the brightness of the galactic 
background 


t(») = f K«)d.e 


The observed brightness at the spacecraft is therefore 

-T(!) 


— t ( co) 

I(o) = I (galactic background) e + 


(a 

•'o n 


JO) 


d£ 


(A) 


(1.15) 


15 



GALACTIC BACKGROUND 



Figure 1.5 AN EXAMPLE OF A RAY PATH IN THE GEOMAGNETIC EQUATORIAL 
PLANE ALONG WHICH EQUATION (1.11) IS TO BE CALCULATED. 


If the brightness due to the galactic background is subtracted from 
the total observed brightness, the contribution from the magnetosphere 
is given by 

r 1 "T(i) T ( . 

I (magnetosphere) = / ^ — d£ (1.16) 

•'o n (jO 

r 

+ 

This equation indicates that the brightness along a ray path through 

the magnetosphere is the sum of the brightness contributions from 

all interior points of the emitting region along the ray path reduced 

by their appropriate factors, e , which account for the reabsorption 

of the radiated energy. The units, steradians, in the definition of 

2 

brightness and the factor n in the preceding equations account, 

+ 
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respectively, for the source region subtending a finite solid angle 

at the point of observation and for the focusing and defocusing of 

the electromagnetic energy as it propagates through the medium. 

To evaluate Equation (1.16) we must determine the emission 

coefficient J(4) and the absorption coefficient K(£) , both of 

which are functions of the energetic electron population in the 

2 2 

magnetosphere, as well as n (n , in turn, depends on the thermal 

r + r + 

plasma) , Once these quantities have been determined along all ray 
paths passing through the emission region and terminating at the 
observer, an average brightness (I ) for the magnetosphere can be 
calculated. From I , the power received in a frequency interval 
Af^ about f^ can be calculated from 

AP r (Af x ) * i *A e (f 1 ).fi(f 1 )^(Af 1 )^I av (f 1 ) (1.17) 


where 

I (f ) = average brightness of the magnetosphere at frequency f 
a v i i 

(defined in Section III,G. 2) 

A (f , ) = effective area of the spacecraft receiving antenna at 
ex 

frequency f^ 

Q(f^) = solid angle subtended by the emission region at the 
satellite at frequency f^ 

In this equation it is assumed that the receiving antenna is isotropic. 
When this is not true, the antenna beam pattern and the variation in 
brightness with observation direction must be considered [Kraus, 1966] . 
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D. Gy ro- Synchrotron Radiation from a Single Electron 

To obtain a representative value for I (magnetosphere) , it is 
necessary to understand the properties of the gyro- synchrotron 
mechanism. To this end, we will develop the more complex case of 
emission in a magnetoplasma by first considering the more fundamental 
situation of radio noise generated by electrons spiraling along 
magnetic field lines in free-space. As a starting point, the model 
for the particle executing circular motion about the field line will 
be presented (i.e., a, the particle pitch angle defined in Figure 
1.3, will be assumed equal to ninety degrees) and then generalized 
to the case of helical motion. Finally, radiation in a cold plasma 
will be discussed as an extension of the simpler free-space results. 

1 - Radiation from a Superthermal Electron in Free-Space with g = 90° 
For a moving charge to produce electromagnetic radiation in 
free-space, it must have a non-zero component of retarded acceleration 
with respect to the observer. When the acceleration has a component 
directed opposite to the particle velocity, the particle is slowed, 
and bremsst rahlung , or braking radiation, is emitted. When the 
charge is moving in a static magnetic field, the velocity and 
acceleration are perpendicular, and the resulting radiation is called 
magnetic bremsstrahlung, or gyro-synchrotron radiation. 

For particle energies E> 100 kev, the radiation is concentrated 
in the orbital plane of the electron. Furthermore, the energy is 
radiated in an increasingly narrow beamwidth as the particle energy 
increases, and the emission is directed in the particle 1 s instantaneous 
direction of motion. For example, the total power radiated by an 
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electron with JE 100 kev has a radiation pattern whose beam area is 

approximately 0.82 steradian, while for E — 1000 kev the value is 

0.10 steradian. Thus, an observer fixed with respect to the static 

magnetic field will detect a sequence of pulses separated in time by 

T = 1/f = y/f (f is the nonrelat ivistic gyrof requency of the 

g T e e 

particle and Y is the Lorentz factor defined as the ratio of the 

particle’s total energy to its rest energy). The length of each 

pulse is given by At = m e /eBY 2 [Panofsky and Philips, 1962]. 

From this discussion, it is evident that the distribution of the 

radiated power in the frequency domain consists of harmonics of the 

angular relativistic gyrof requency (f^) . This observation can be 

generalized by expressing the power radiated by a single particle, at 

a frequency equal to the s th harmonic of its gyrof requency , as a 

function of five parameters. This function is written as 

P (f E,R,0,n ) watts ster 1 where the variables have the following 
s + 

interpret at ions: 

f, the frequency of observation, where f = sf g 
E, the energy of the particle 

R, the instantaneous radius of curvature of the particle trajec- 
tory. (R is itself a function of the magnetic flux density (B) , 
the velocity of the particle, and the angle (a) which the veloc- 
ity vector makes with the local magnetic field.) 

0, the angle between the wave propagation vector and the 
static magnetic field 

n +f the real parts of the complex refractive indices of the 
medium 


19 



Before discussing the effects of these parameters on P it is 

necessary to obtain the equation for the total power radiated by an 

electron via the gyro- synchrotron process when a = 90° and n =1*0. 

+ 

This problem was considered by Panofsky and Philips [1962], who derived 
the following expression: 


2 2 2 2 r- 

e u) B s 


P s (f,E,R,0) = ^ 


8tt e c 

O 


j' 2 (sp sin 0) + COt / 9 ? - J 2 (sp sin 0) 

$ 3 

a. is) 


where 

e = electronic charge 

m = rest mass of electron 
e 

- permeability of free-space 

e o = permittivity of free-space 

c = velocity of light 

v = velocity of particle 

B = magnetic flux density 

J = Bessel function of the 
s 

first kind of order s 

J 1 = derivative of J 
s s 

f = 00/217 = frequency 
fg “ 00g/2n = relativistic 
gyrof requency 


H = magnetic field intensity 

2 i 

P = v/c, y - l/(l-p ) 

R = radius of circular orbit = 

m vy/eB 
e 

U)- = angular gyrof requency = 

& 

v/R = eB/m Y 
e 

Kj = total energy of particle = 

2 2 
me = ym c 
e 

2 

E = kinetic energy = m e c (y - 1) 

_ 2 
\ - ■„<= 

uo = angular wave frequency 
s - harmonic number 
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Using the above definitions, Equation (1,18) can be recast as 


4 2 2 
e P B 


p s (f,E,B,e> = 2 2 2 


8n Y m e c 
e o 


s 2 j'^(sp sin 9) + 5 C ^ -— ^ ’J 2 (sp sin 9) 

s P 


(1.19) 


(i) To investigate the effect of the particle energy, E, on the 

emission spectrum, P (f,E,B = constant, 0 = 90°) is plotted as a 

s 

function of s for different values of E in Figure 1.6. This 
figure shows that, as E increases, the harmonic number (and there- 
fore frequency) at which peak emission occurs (s ^ also increases. 
In the energy range considered in this study ( E> 100 kev) , s 

pK 

3 

is given approximately by s = [(3/2)y ] , where [x] denotes the 

pk 

largest integer less than or equal to x. 



Figure 1.6 ENVELOPES OF P s (f,E,B,0) VERSUS HARMONIC NUMBER FOR ELEC- 
TRONS OF VARIOUS ENERGIES WHEN Q = 90° AND B = 5.36 x lfl' 
GAUSS (i.e., f e = 15 kHz). 
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(ii) The effect of changes in the magnetic flux density is easily 

seen from Equation (1.19), which shows that P (f,E = constant; B, 0 = 

s 

o 2 

90 ) « B . 

(iii) The radiation pattern for the particle emission is given 
by the polar plots for P g (f ,E = 0.1, 0.5, 1.0 Mev; B = constant, 0) 
in Figures 1.7 a-c. A number of important results should be noticed 
from these curves. (The emission is taking place in free- space; 
therefore, all harmonics are generated and propagated. To keep the 
figures legible, however, all of them were not drawn.) 

(a) The power received for s > 3 and E > 100 kev is a 
maximum for 0 = ce = 90°. 

(b) The beamwidth (defined as the angle between half-power 
points measured about the direction of maximum emission) decreases 
as harmonic number increases . 

(c) Higher harmonics contribute very little power for 
directions near or along the magnetic field line. Lower harmonics, 
however, contribute appreciable energy close to the field lines. 

In particular, the first harmonic is nearly twice as strong along 
the field line as compared to 0 = 90° (see the E ^ 100 kev plots) . 

(d) There is symmetry with respect to 0 = 90°. (This fact 

is readily deduced from Equation (1.19) # ) - 

(e) The power received at a frequency f^ = sf^ is strongly 
dependent on the angle of observation 0 . 
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POWER twotts/sterodion) 


B 


Figure 1.7 a-c P s (f,E,B,8) VERSUS THETA 

FOR VARIOUS HARMONIC NUMBERS AND 
THE VALUES OF THE PARAMETERS 
GIVEN IN EACH FIGURE. 



Figures 1.7 a-c and Equation (1.19) represent the total power 
radiated by the electron, which is defined as the sum of the power 
radiated into the ordinary and extraordinary modes. These modes of 
propagation will be defined and discussed in Section l.D, 4. 

2. Radiation from a Superthermal Electron in Free-Space with a ^ 90° 
The model presented thus far is the simplest case for which 
gyro- synchrotron emission can occur, and it has been discussed and 
applied extensively in the literature [Peterson and Hower, 1966; 

Hower and Peterson, 1964; Vesecky , 1967]. However, if the particle 
travels in a helix along the magnetic field line (i.e., a £ 90° with 
n+ = 1.0), the model becomes considerably more complex. 

For this case, the appropriate equation can be derived by 
integrating Equation (6.14) of Bekefi’s [1966] with respect to radian 
frequency (w) . The result for the total radiated power per unit 
solid angle is 

P s (f,E,B,0,a) = 

L' 2 / 5 ?! 5 ™ 6 V /c os 0-p 2 \ 2 j2/s(3 lS in 6 \ 
8n 2 e o cmJ Y 2 (l - p 2 cos Q) 3 “l 1 ^ 0080 / \Pi sin Q / * \ l-0 2 cose ) 

( 1 . 20 ) 

where the frequency of observation (f) is now given by 

f = S f r i i a. 2 i) 

g i - p 2 cos e 


24 



and 


P 1 = P sin o' 
P 2 = P cos or 


( 1 . 22 ) 


The term (l-p cos 9) 1 is due to the Doppler shift occurring as a 

A 

result of the electron motion with respect to the fixed observer. 

Calculations by this author show that the effects of E and B 

on P (f,E,B,0,cv) satisfy sub-paragraphs b,c,d and e of the preceding 

section where now f = sf is replaced by Equation (1.21), For the 

g 

present case, however, the maximum radiation for any harmonic and 
energy is in the direction of the instantaneous orbit of the electron 
(as shown by Figures 1.8 a-c for a = 70 and Figures 1.9 a-e for 
a = 30°). Furthermore, the symmetry about 0 = a no longer exists, as 
the power falls off much more rapidly for 9 < O' than for 9 > or . 

This effect is a result of the modes of propagation alluded to in 
the last paragraph of the preceding section. As will be discussed 
later, the majority of the energy in the "ordinary” mode is generated 
near 0 = 90° while the energy in the "extraordinary" mode is generated 
predominantly along the direction 0 = a . The sum of the power 
radiated into both of these modes, therefore, yields the skewed 
plots shown in these figures. 

The curves also indicate that the peak power for a given harmonic 

o 

decreases as the pitch angle decreases from 90 . This phenomenon is 
clearly depicted in Figure 1.10, where each curve is plotted for 
0 = a. Over the energy range 100 kev < E < 3 Mev, calculations showed 
that the peak power for a given harmonic was within a factor of one 
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to 

cr. 



Figure 1 . 8 a-c P (f s E,B,9,cr) VERSUS THETA 
s 

FOR VARIOUS HARMONIC NUMBERS AND 
THE VALUES OF THE PARAMETERS GIVEN 
IN EACH FIGURE.. 


FREE-SPACE 100° 80° 




to 


Figure 1.9 a-c P (f ,E,B,6,a) VERSUS THETA 
s 

FOR VARIOUS HARMONIC NUMBERS AND 
THE VALUES OF THE PARAMETERS GIVEN 
IN EACH FIGURE. 




half from that at cr = 90°, when 
than 70°, the peak power begins 


o o 

90 < a < 70 . 

to decrease at 


As a becomes less 
an increasing rate. 



Figure 1.10 ENVEI/DPES OF P (fjE.B^.o) VERSUS HARMONIC NUMBER. EACH 

s 

CURVE IS OBTAINED FOR THE OBSERVATION DIRECTION COINCIDING 
WITH THE INSTANTANEOUS ORBITAL PLANE OF THE ELECTRON 
(i.e. , 6 = a ) . 
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This section has extended the gyro-synchrotron model to include 
particles whose pitch angles are not equal to ninety degrees, which 
is the general case applicable to the free-space problem. It remains, 
however, to investigate how the power spectrum changes when the 
particle radiates in a medium where the refractive index is no longer 
equal to unity. 

3. Radiation from a Superthermal Electron in an Isotropic, Refractive 
Medium 

Before discussing the final case, we shall consider an isotropic 

medium composed of a tenuous plasma satisfying the conditions derived 

in Section I.B (i.e., the magnetic field is assumed weak enough not 

to affect the propagation of electromagnetic waves yet strong enough 

to allow gyro-synchrotron radiation to occur). A medium of this 

nature is described by one refractive index, n , for both the ordinary 

and extraordinary modes of propagation. Under this condition, the 

* 

power spectrum P (f , E,B ,0 ,a ,n ) can be obtained by modifying the 
s 

derivation for n + = 1.0 in Bekefi [1966], This modification is based 
on the fact that the phase velocity of the electromagnetic waves 
propagating in the plasma is not c (the velocity of light in free- 
space) but is instead 


v 

ph 



(1.23) 


Using this equation in Bekefi 1 s analysis yields 
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* e s B p n 

P s( f ,E 'B>9 >Q'> n ) = 2 2 2 * 3" 


8n e cm y (1-n 8 0 cos 6) 
o e r & 


3^£ 

2 / sn P 1 sin e \ 

S yl-n^p^cos 0 / 


^ 2 ^ 

( cos e-n p \ 2 / sn p^sin 9 \ 

n^sin 0 / S \ l-n*p 2 cos 0/ 


(1.24) 


where f is now given by 


f 



1-n 


* 

P 


1 


cos 0 

A 


(1.25) 


This formula was evaluated for n =0,7 and ot — 90°, and the results 
were plotted in Figures l.lla-c, These figures do not represent a 
realistic case as n is actually a function of frequency (and thus 
harmonic number). They do, however, emphasize the relationship 
between p s ( f > E » B ,0 9 ct ,n ) and n . As seen, for increasing particle 
energies, the power radiated by a particle is severely decreased in 
comparison with the n = 1.0 case. This effect is so pronounced that, 
on a scale consistent with previous figures , only the tenth harmonic 
for E s 1.0 MeV can be drawn. These figures also show this phenomenon 
to be greater at higher harmonics , as evidenced by comparing Figures 
l.lla-c with Figure 1.7a-c. These same observations are also applicable 
to pitch angles other than 90° when n* is less than unity. 
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* 

Figure 1.11 a~c P g (f ,E , B,0 , a,n ) VERSUS THETA 
FOR AN ELECTRON RADIATING IN AN 
ISOTROPIC MEDIUM WHERE IT WAS 
ASSUMED n (f) =0.7. 



Even though the above analysis helps to clarify the effects of 
* 

varying n and P , it does not represent a physical situation. To 
s ‘ 

understand what actually occurs, consider the emission from a super- 

thermal electron in a magnetoplasma with a plasma frequency (f ) of 

P 

15 kHz. Since the medium is assumed isotropic (i.e., the effects of 
the magnetic field are ignored) the refractive index is given by 
[Tanenbaum , 1967] 



Using this equation for n in Formula (1,24), we calculated the curves 

in Figures 1.12a-c, which show a number of new features. 

First, both the power radiated at the lower harmonics and the 

beamwidth of these harmonics decrease in comparison to the free-space 

problem. However, as the harmonic number increases, these effects 

become less evident. This phenomenon was first noticed by Razin [i960] 

and is easily explained by noting that as the frequency (and hence 

* 

harmonic number) increases, n approaches unity, its value for free- 

space. Razin has shown that the medium may be treated as free-space 

provided the frequency of the electromagnetic radiation (f ) satisfies 

2 

the relationship f » 3/2 f /f 

P e 

Second, a number of the lower harmonics for each energy are no 
longer emitted. (See Figures 1.12a-c, where the lowest harmonic 
plotted is the first to be radiated.) As shown, the greater the 
particle energy the larger the harmonic number which is first 
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Figure 1,12 a-c (f ,E ,B,0 , a ,n ) VERSUS 

THETA FOR AN ELECTRON RADIATING 

IN AN ISOTROPIC MEDIUM WHERE n* 

IS A FUNCTION OF FREQUENCY AND 

f = 15 kHz. 

P 


n 100° 80° 



generated. This phenomenon is explained by noting that for emission 

* o 

and propagation to occur, n must be greater than zero. For a = 90 

. * 

and f = sf “ sf /Y , n can be expressed as 
g ^ 


* 


n 



(1.27) 


Thus , 



(1.28) 


This expression is plotted in Figure 1.13, where the particle energy 
is given as a parameter. This plot only has values at integral values 
of s ; the solid curves are drawn to indicate the envelope of 
Inequality (1.28). 

For pitch angles other than ninety degrees, the features described 

in the last paragraph still apply, as seen by comparing Figures 1.14a-c 

o * 

with 1.8a-c. For ^ 90 , however, n is an implicit function of 0 
and . That is, by substituting (1.25) into (1.26), n* is now given by 


* 

n = 


- f cos e + sf 

p 2 g 


J 


~ 2 a 2 2 2 2 2 
f cos e + s f - f 

p 2 g p 


2.2 22 2 
s f + f B_ cos 9 
g 


(1.29) 
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Figure 1.13 ENVEIXJPE OF THE REFRACTIVE INDEX, n* , VERSUS HARMONIC 

NUMBER FOR VARIOUS ELECTRON ENERGIES. THE MEDIUM IS 

ISOTROPIC AND IS REPRESENTED BY f =15 kHz and 

e 

f =15 kHz* 

P 

This result allows for the interesting possibility that a given set 
of parameters (s,f^,(3 2 and f ) will result in n becoming less than 
zero for small values of the angle theta. The consequence of this 
situation will be demonstrated in the following section where the 
theory is extended to the general case. 
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4 . Radiation from a Superthermal Electron in a Cold, Anisotropic Plasma 


Radiation from a superthermal electron in a cold, collisionless, 
magnetized plasma has been investigated by several authors [Eidman, 1958 
McKenzie, 1967; Melrose, 1968; Truelson and Fejer, 1970] who applied 
their models to explain VLF and LF radio noise observed at ground- 
based stations. As a result of their interest in low frequency pheno- 
mena, these authors were primarily concerned with whistler mode 
propagation and hence regions in the magnetosphere where the refractive 
indices, n + , (at VLF and LF frequencies) were substantially greater 

than unity (although their results were valid for n ^ 0) . The 

+ 

model that gives the emission in an arbitrary anisotropic, collision- 
less, cold medium is therefore obtained from these studies, but is 
recast in terms of the parameters X and Y defined in Figure 1.2. The 
pertinent equations are [Melrose, 1968] 


P (f ,E,B,0 ,<*,n ) = 
s + 


4 2 2 2 

e s B n.T. 

JJ 


2 2 2 2 
8tt e o cm e Y (1 - n jP 2 cos 6) 


l l J s <5) -( 




k«> 


1 - B„n .cos 9 
2 J 


w dn . \ 

1 + ITT — ji ) 

J do) / 


(1.30) 


in watts/steradian, where 


f = sf 


1 - n.p cos 0 
J & 


(1.31) 
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(1.32) 


! 


sn R sin 0 
3 ^ 

i “ n jp 2 cos 6 


A - H.cos 0 + T sin 0 (1.33) 

y 3 A 


A = T .cos 0 - H.sin 0 
z J J 



1 + H 



(1.34) 


(1.35) 


h 

1H. = ^ * 1 

J h 


2Y(l-X)cos 0 


2 2 — f 4 4 2 - 2 2 *" 

Y sin 0 +s/Y sin 0 + 4Y (1-X) cos G 


(1.36) 


E Y sin 0 + H.Y Xcos 0 sin 0 

iT = — — i - 

j E ,.2 i . . . „2 2 . 


(1.37) 


-x 1 - Y“ - X(1 - Y cos 0 ) 


n . = 1 - 


2X(1-X) 


2 (1-X) - y 2 sin 2 0 + [y 4 sin 4 0 + 4Y 2 (1-X 2 ) cos0] /2 


(1.38) 


If no restrictions are placed on the magnitudes of n + (other 
than n + > 0) there is the possibility that anomalous gyro- synchrotron 
(given by s < 0) and Cerenkov (given by s = 0 ) radiation can be 
generated concurrently with normal gyro- synchrotron emission (s > 0) . 
Result three of Section I.B, however, stated that for the regions of 
the magnetosphere which we shall be investigating, the frequency of 
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the radio noise which is generated will be above the local upper 

hybrid frequency (f u ) . This fact implies that f will also be 

greater than the ambient electron gyro and plasma frequencies. 

Therefore n £ 1 which, when taken together with Equation (1.31), 
* + 

gives 


cos 6 



Vs 


<2 1 


(1.39) 


Solving for sf^/f 

sf 

—Iso (1.40) 

f 

and 

s s 0 (1.41) 

However, for s = 0, sec 0 = (n p ) ^ 1, and hence n S (0 ) . In 

general, p <1 resulting in n necessarily being greater than unity 

for s = 0. The medium under study will therefore give rise only to 

normal gyro- synchrotron emission (s > 0). In this section, this form 

of radiation from energetic electrons immersed in the complex medium 

described above will be investigated. 

Returning to Equation (1.30), the total power radiated at a given 

harmonic is the sum of that emitted into two modes (represented by 

the summation over j), which are the ordinary (n j = 1) and the 

+ 

extraordinary (n , j = 2) modes defined by Ratcliffe [1962]. To 
discuss the total power in general terms is intractable due to the 
complexity of the mathematics. Nevertheless, a few observations can 
be made by noting that the general phenomenon is only an extension of 
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the simpler cases discussed previously. Calculations made by this 

author showed that the basic characteristics listed in the preceding 

* 

sections (i.e., the effects on P (f , E,B,0 .o^n ) of changes in its 
parameters) remain basically the same for each mode when they are 
considered separately. Thus the total flux, being the sum of the 
power radiated into each mode, demonstrates the same dependence on 
these parameters . 

More specifically, Equation (1.30) was used to calculate the 
emission expected from a superthermal electron under geophysical 
conditions and locations in the magnetosphere pertinent to this paper. 
The parameters for two such cases are listed in the headings of 
Figures 1.15a-h and 1.16a-h respectively. The reason for this choice 
of parameters and how they relate to this study will become evident 
in the next chapter. It is intended at this point only to show the 
effects of the complex medium on the radiation. The new results are 
summarized below: 

a. The maximum power emitted in the ordinary mode is a factor 
of two or more below that generated in the extraordinary mode for all 
electron energies E. This ratio is a function of E, increasing as 
E increases. (In free-space, however, calculations for E ^ 10 kev 
showed that this ratio lies between 0.3 and 3.0 where 1 0 | s (o' - 10°| 
and increases rapidly for [e| = \a - 10° | . Thus, for n + = 1.0 and 
for particles whose instantaneous orbital planes are directed near 
or along the direction of observation, the power radiated into the 
ordinary mode is comparable to that radiated into the extraordinary 
mode . ) 
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Figure 1.15 (continued) 





5<Kp<6" 
ORDINARY MODE 
£=IOO kev 
a =50° 

f e =20kHz , 

fp=20 kHz X 

i4oy 


10 0° 80° 


ieo Q i^ i 
I0' 27 
(e) 



IO' 30 fO" 30 I0" 29 

POWER (wotts/steradion) 


5<Kp<6 
ORDINARY MODE 
E = I Mev 
a = 50° 
f e = 20 kHz 

f n = 20 kHz / 


■ > i ti i i 

I0“ 27 I0' 28 


100® 80° 



10' 29 10" 30 IO" 30 

POWER (watts/steradian 


5<Kp<6 + 

EXTRAORDINARY MODE 
E = 100 kev 

a = 50° ^OV" 

f e =20kHz X"' 

fp=20 kHz X 
140 Y 



IQ-50 j0 -5O | 0 -Z9 

POWER (wolls/slerodion) 



5<Kp<6 + 

EXTRAORDINARY MODE 
E- I Mev 

a = 50° I20X 

f e =20 kHz 
fp-20 kHz X 
I40X 



I O ' 50 I O ' 30 IO ' 29 

POWER (watts/steradian) 


Figure 1.16 (continued) 


b. The maximum power in the extraordinary mode is radiated 


predominantly along Q = or, while the ordinary mode is emitted closer 
to 0 = 90°. This feature is particularly evident for the lower energy 
electrons and lower harmonic numbers . If these two modes are summed 
to obtain the total power radiated into a given harmonic, the resul- 
tant polar plot for P (f ,E ,B ,9 n ) will be asymmetrical, as described 

S 4- 

in Section I.D.2. This result is general and applies to the free-space 

problem since Equation (1.30) reduces to (1,20) when the appropriate 

free-space values for n ,T. and H. are used in (1.30) [Eidman, 1958] 4 

J J J 

c. The Razin effect (Section I*D.3) is evident for each mode, 
as seen by comparing these figures with those of previous sections. 

d. The lowest harmonic generated (i.e., n > 0) is the first 

+ 

plotted in each figure. Furthermore, the lower harmonics can be 
abruptly cut off for small values of the angle theta. This is a 
result of Equation (1.31) not being satisfied for values of theta 
less than the critical value 0^ given by 


e 

c 



(1.42) 


in which the following inequality must be satisfied; 


sf 


1 - 


<; i 


(1.43) 


n J 


+ K 2 
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That is, for a given set of parameters s, f , f and values of 0 

less than 0 result in values of n (n being a function of f . f. f 
c + + g p 

and 0) which do not allow Inequality (1„43) to be satisfied. 

These results conclude the presentation of the radio noise model 
which shall be used in the subsequent text to quantitatively determine 
and discuss the noise generated in the magnetosphere. To this point, 
however, we have only determined the model for the power radiated at 
the site of its source, a single superthermal electron. We must now 
turn to the problem of theoretically deriving the total energy 
radiated by an ensemble of energetic electrons and the subsequent 
propagation of this energy from the source region to the observer. 

E. The Power Spectral Density Emitted from a Single Electron 

In the previous sections the characteristics of P (f ,E ,0 ,<* ,n ) 

s + 

as a function of its six parameters were investigated in detail to 

give the reader an understanding of the basic properties of the gyro- 

synchrotron mechanism. However, to calculate the emission coefficient 

J(4) it is necessary to derive an expression for the power spectral 

density P(f , E,B,0 ,a ,n + ) which is expressed in units Of watts Hz 1 
-1 

steradian 

The derivation of P(f ,E,B,0 ,<y ,n ) as a function of the six 
variables proves to be unnecessary. As stated in paragraph eight of 
Section I.B (and justified in Chapter II) , the magnetosphere shall 
be modeled as a stratified medium composed of concentric shells. 

Within each shell at the geomagnetic equator, B and can be taken 
as constant. Furthermore, once the position of an observer relative 
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to the earth is specified, the parameter 0 and hence a can also be 
taken as constant within each of these shells. (The angle <x is 
determined since emission in the direction 0 results from particles 
with a & 0 due to the directivity of the gyro- synchrotron mechanism.) 
Thus the general function P(f ,E,B,0 ,a reduces to evaluating 
P(f,E) for given sets of values for 3,0,0, and n . 

Theoretically, P(f,E) is the power radiated by an electron, with 
energy dE about E, into a bandwidth df about f. In Section I*D.4, 
we found that for an anisotropic medium, f is given by 


sf 

f = 2 

■y(l - n + p 2 cos 0) 

where 

E + E 
o 

Y “ E 

o 

and by definition 


(1.44) 


(1.45) 


(3 = p cos <y = 1 1 - cos a 

V Y 


(1.46) 


Substituting (1.45) and (1.46) into (1.44) yields 


sf 


E +E 
o 


( 


2 f 2 

s i 

e 2 2 2 

e i — ?r- + n cos 8 cos a 

f~ ± n + cos a cos 0 \ f 2 + 


i 

■ 


(1.47) 


2 2 2 
1 - n cos <y cos 0 
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The + sign in front of the radical indicates that electrons with 
two different energies, E and E can generate electromagnetic 

i 4 

waves at a frequency f, corresponding to a given harmonic, s. 

Differentiating the above expression with respect to f and 
solving for |df| gives 


| df - 


2 2 2 2 
f (1 - n cos & cos 0) 
dfc 


sf + 
e — 


2 2 

s f n.cos or cos 0 
e + 


“ IdE 


, 2 2 2 2 2 2 2'/i 

is f + f n cos 0 cos 0 - f ) 
e + 


( 1 . 48 ) 


Then the power spectral density for a single electron can be 
written as 


P(f ,E) 


p s (f,E) 
[df | 



(y - n_/y 2 -l cos 9 cos Q/)f 


f 

e 


( 1 . 49 ) 


where 

P(f,E) = power radiated into a frequency interval df about f 

by an electron with energy dE about E (watts" Hz 1 ster 1 ) 

P g (f ,E) = power radiated by an electron into its s th harmonic 
-1 

(watts ster ) 

f = frequency of observation 

E = center of energy range within which the electron being 
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considered must be found. Note that this energy is given 


explicitly by Equation (1.44) once f,n ,0 ,f , 0 - and s are 

■f* e 

specified 

6( ) = Dirac delta function which implicitly gives the energy 

E that the electron must have to radiate at a frequency f for 
a given s 

If we consider the specific case where a = 90°, y » 1 and 
s» 1, Equation (1.48) reduces to 


E sf 

|dE| = - — 2 6 | CUE | (1.50) 

f 

From the argument of the delta function 

YE 

|dE| M — (1.51) 

Furthermore, for these values, P^(f ,E) may be considered continuous 
and s ssyf/f^ [Peterson and Hower, 1966; Vesecky, 1967]. Substituting 
these results into (1.49), we obtain 

P _<f ,E) 

P«,E) = (1.52) 

S 

This equation is just the approximation used by previous authors 
[Peterson and Hower, 1966; Vesecky, 1967] when they considered prob- 
lems related to the detection of synchrotron radiation at frequencies 
where s» 1 from relativistic electrons. In this study, however, 
electrons with energies E > 50 kev and s ^ 1 must be considered; there- 
fore, the general Equation (1.49) will be used. 
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F . The Gyro-Synchrotron Emission Coefficient from an Ensemble of 


Superthermal Electrons 

For the problem considered in this paper, the emission coefficient, 
J(4), can be defined as [Vesecky, 1967] 


J(f,e,B,n + ) = JJ P(f ,B,B,e,a,n + ) N(E,<y)dE da 


(1.53) 


where 

J(f,9,B,n ) = power radiated at a frequency (f) and angle (0) 
from the static magnetic field by an ensemble of electrons with 

HHffownt oruar and pitch angles, specified in watts 

-3 -1 ^ -1 

m Hz ster 

P(f ,E,B,0 = power spectral density at a frequency f 
radiated by an electron with energy dE about E and pitch 
angle a at an angle 9 from the magnetic field line , given 
in units of watts Hz *ster *" 

N(E,qO = the number density of electrons with energy dE about 

E having pitch angle do about o, specified in units of electrons 

~3 -1 

m ster 

As stated, the magnetosphere is modeled as a stratified medium in which 

9 ,B, and are taken to be constant within each layer. Furthermore, 

the pitch angle distribution of the superthermal electrons satisfies 
.2 

a sin <y law (see Section II. B. 5) which is much broader in angular 
extent than the synchrotron emission patterns calculated in the 
preceding sections. Consequently, the integration with respect to 
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O' in Equation (1.53) reduces to determining the number of electrons 
which radiate in the direction of 0. Thus, we may express Equation 

(1.53) as 


J(f) = 


/ 


P s (f,E) 
I df I 


N(E) 


s - 


(Y-n + J^ 


1 cos 0 cos a J f 


dE 


(1 ,54) 


where 


N(E) = 2tt 





, N(E r of) sin a da 


(1.55) 


In the last equation, the limits of integration and a 2 are a 
function of 0. They will be discussed in depth in Sections III.G.3 
and III.H.2. 

By virtue of the delta function, Equation (1.54) reduces to the 
following summation 

V V f ’ E > 

j(f) « 2, a 

s=l | df j 

where the delta function is retained in this and subsequent equations 

to show that the particle energy (E) is determined once a set of 

values for s,n ,f ,0 t <y and f is specified. Equation (1.56) gives 
4" e 

-3 -1 -1 

the power in watts m Hz ster radiated at a frequency (f) from 


s - 


{y~n +f Jy "1 cos 0 cos * ) 1 


(1.56) 
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an ensemble of electrons located in a layer of the magnetosphere 
represented by the quantities B,0 and n . 

This equation, however, is not in the most convenient form for 

this study, as the particle density data, given in the next chapter, 

-2 -1 -1 

is expressed as a flux (electrons m kev sec ) . It is therefore 
more useful to write (1.56) in the following form 

| y-n+Jy 2 -l cos 0 cos <y ) f 

s — 

f 

e 

(1.57) 

where |dE/df] is given by Equation (1.48) in kev Hz 1 and the particle 
flux N f (E) is given by N p (E) = c(3(E)N(E) per kev. 

G. Absorption Mechanisms and the Absorption Coefficient, K(l) 

As the electromagnetic waves propagate through the magnetosphere, 
it is possible that they can lose energy to the medium from collisions 
between electrons and other particles [Ratcliffe, 1962] or from gyro- 
synchrotron self -absorption [Pacholczyk, 1970], In paragraph five 
of Section l.B, it was shown that the first process may be neglected. 
The second, however, must be considered in more detail. 

We proceed by using the result obtained by Pacholczyk, who gives 
the absorption coefficient at a frequency (f) as 


J(f) 


-2 

S=1 


N (E) 

if 1 4(E) i 
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dE 


(1.58) 


K<f> = 7 ■[ ^a[ A H 

where 

K(f) = power absorbed by electrons from an incident electro- 

A 

magnetic plane wave propagating in a direction I with a given 
polarization P (nepers m *) 

A 

A(f ,E) = energy radiated at a frequency (f) in the direction I 
with a given polarization P by electrons with energy dE 
about E (watts Hz * ster *) 

This equation can be integrated by parts to yield 


K(f) = — ) N(E) A(f ,E) 


A(f 5 E) 


(1.59) 


The first term in this equation vanishes since N(E) — >0 as E —Voo, 
and A(f,E) — ^0 as E->0. Furthermore, from the definition of both 
terms, A(f,E)= P(f,E) and 



Substituting (1.49) into this equation and once again retaining the 
Dirac delta function, we obtain 
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K(f) = - 


s 


C 2 y 2 P s (f ' E) d N(E) ~| 

f 2 sTi |d, | dE [_ E 2 _ 


6 


(Y ~ Jy • -1 cos 6 c °s cy)f 
f 

e 


( 1 . 61 ) 


This last equation and the one obtained for J(f) will be used in 
Chapter III to obtain numerical solutions to the energy transfer 
equation defined in Section I.C. Before this is possible, however, 
appropriate models for the magnetosphere must first be developed. 
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CHAPTER II 


MAGNETOS PHERIC MODELS AND RAY PATH CALCULATIONS 
A . Introduction 

In order to calculate the gyro- synchrotron radiation generated 
in the earth’s outer radiation belt and propagated to an observer in 
the interplanetary medium, we must develop models for the energetic 

and thermal electron populations in this portion of the magnetosphere, 

\ 

This problem is complicated, however, by the dynamic behavior of these 

quantities resulting from geomagnetic activity. An attempt will 

be made, therefore, to develop models that are representative of 

average conditions for these populations during three levels of 

geomagnetic activity: 1 ^ ^ 2, 3£K^^4, and 5 ^ K p ^ 6 , 

where K is the planetary magnetic three- hour- range index. These 
P 

models will then be used with the radio model presented in the 
previous chapter to calculate the intensity, spectrum and variation 
of the noise generated in the magnetosphere as a function of geomag- 
netic activity. Finally, these theoretical results for the noise 
will be compared to the data returned by the GSFC radio experiment 
on the IMP-6 spacecraft. 

These particle models are to be developed from twp sources of 
data: open literature and personal communications. One difficulty 

arises, however, since the data obtained from these sources were 
collected prior to the time period covered by the results from the 
IMP-6 experiment. These data, pertaining to the particle populations, 
will therefore be extrapolated in time to allow a direct comparison 
between the theoretically calculated and experimentally measured 
radio noise. 
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B. The Superthermal Electron Population 

The change in the density of trapped electrons in the outer radia- 
tion belt due to geomagnetic activity is believed to be caused by an 
acceleration and subsequent injection of low energy electrons into 
the stable trapping regions [Williams and Smith, 1965] . The source 
of these low energy electrons is not entirely clear at present. 

However, it is found that the effect of a sudden intensive magnetic 
disturbance is an increase in trapped electron densities for L> 3. 

(In this study, Mcllwain's magnetic coordinates [Mcllwain, 1966] will 
be used to discuss the location of trapped electron fluxes in the 
earth 1 s magnetosphere.) 

Frank et al [1964] found this impulsive increase in the elec- 
tron intensities to follow an orderly sequence. During the main 
phase of a magnetic storm, the density of low energy electrons, 

40 kev £ E £ 200 kev, begins to increase within one day after the 
onset of the storm. Following a delay of several days, the intensity 
of electrons with E ^ 1.6 Mev also increases over pre-storm values 
by a factor similar to the enhancement of the low energy component 
of the fluxes. For these fluxes in the outer radiation belt this 
factor can be 1 to 1.5 orders of magnitude over quiet time flux 

intensities for K & 3 (moderate activity), and as much as 1.8 to 
P 

2.0 orders of magnitude for « 5 [Freeman, 1964; Frank et al , 1964; 

Owens and Frank, 1968; Craven, 1966; Williams, 1966; Rose, 1966]. 

Frank further observed that the intensity of electrons, 40 kev ^ E 
£ 200 kev, decays over a period of a week from its maximum, disturbed 
value to its pre- storm level; the more energetic electrons take a 
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considerably longer period of time. The temporal history for the 
intensity of electrons 200 kev < E < 700 kev is intermediate to the 
low and high energy components of the fluxes. 

These large, temporal variations in the trapped electron popula- 
tion make it difficult to determine an average value for the terres- 
trial noise that is generated by these particles via the gyro-synchro- 
tron mechanism. To simplify the problem (and not neglect these 
particle fluctuations), models for the trapped electrons as a function 
of will be developed. That is, we shall determine the expected 
increase in the electron densities during periods of sustained 
geomagnetic activity. As pointed out in the introduction, these 
models are arrived at by using two sources of data, results on the 
trapped particle population published in the literature and data 
obtained through personal communications with Dr. Harry West (at 
Lawrence Livermore Laboratories, Livermore, California) who supplied 
published and unpublished results from the OGO-5 electron spectrometer 
experiment. 

1 . The OGO-5 Electron Spectrometer 

A very brief description of the electron spectrometer aboard 
the OGO-5 is presented to familiarize the reader with this experiment. 
For a more detailed discussion of this and other experiments on the 
spacecraft, see West et al [1969]. 

The electron spectrometer is composed of seven discrete channels 
that respond to electrons whose energies lie in the ranges given in 
Table 1. This table also includes the acceptance angle for each 
channel, defined as the effective width at 50% acceptance in the 
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scan plane [West et al , 1972]. Pitch angle data for the electrons 

were also obtained in this experiment by scanning the entire aperture 
of the spectrometer [West et al, 1969], 


Table 1 

DESCRIPTION OF THE CHANNELS COMPRISING THE ELECTRON 
SPECTROMETER ABOARD THE OGO-5 SPACECRAFT 


Channel No. 

Center Energy in kev 

Acceptance Angle 

1 

79 

7.6° 

2 

158 

5.9° 

3 

266 

4 , 7° 

4 

479 

3.5° 

5 

822 

5.3° 

6 

1530 

4.1° 

7 

2820 

2.5° 


By combining the data received from this experiment with that 
available from the literature, it was possible to obtain a reasonable 
estimate for the differential energy spectrums for the electrons 
trapped in the outer radiation belt. This model will be assumed 
representative for the period of time covered by the IMP-6 radio 
noise data. 

2 . Differential Electron Flux Spectrums for 1 ^ K p ^ 2 in the 

Geomagnetic Equatorial Plane 

Vette et al [1966] have developed a trapped electron environment 
for both the inner and outer radiation zones through L = 6.0, Their 
model is an average representation of the energetic electron population 
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for the time around August 1964, a period of little average geomagnetic 
activity during solar cycle minimum. This model, termed the AE2 
electron environment, was used to generate the differential electron 
flux spectrums shown in Figure 2.1. They show the average perpendicular 
electron fluxes that were measured or extrapolated to the geomagnetic 
equatorial plane. (Vette's model was expressed in terms of omni- 
directional fluxes, which were converted by this author to perpendi- 

2 

cular fluxes to allow comparison with OGO-5 data. A sin CC pitch 
angle distribution was assumed in the conversion. This assumption 
will be justified in Section II, B, 5. 

However, there is one problem with the AE2 model; for L > 5 at 
the geomagnetic equator, in situ data sufficient for an accurate 
representation of the fluxes could not be obtained. As stated by 
Vette, "extrapolations to higher L values should be interpreted with 
caution and only used for order of magnitude estimates” [Vette et al , 
1966], Since the region 5 £ L £ 7 is important to this study, Vette's 
data will be supplemented with in situ measurements made by OGO-5 
for 2 £ L £ 9 and < 1 (Figures 2,2a-b). A comparison of Figures 
2.1 and 2.2b shows that they agree very well for L <. 5.0 and 
E £ 500 kev. But, for L £ 5.0 and E> 1.0 Mev, the OGO-5 data give 
fluxes that are larger than Vette's results by factors as great as 
eight. This difference may be due to the solar cycle; the OGO-5 
data were taken during a time of maximum activity (1968). However, 

Vette [1971] states that as the solar cycle goes from minimum to 
maximum, an average increase by a factor of four is expected for 
electron fluxes with E « 500 kev. Since the fluxes at E « 500 kev 
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are the least sensitive to geomagnetic activity [Frank et al , 1964; 

Freeman, 1964], both the lower, 40 kev £ E £ 200 kev, and the higher 

energy electron densities, E> 1,0 Mev, should show increases larger 

than this factor. Contrary to this expectation, Figures 2.1 and 2.2b 

show that the 0G0-5 and AE2 data agree almost exactly for all L values 

and E £ 500 kev. The implication is that this set of OGO-5 measurements 

shows fluxes lower than those expected during solar cycle maximum; 

therefore, the discrepancy between these Figures (2.1 and 2.2b) for 

h > 5 and E £1.0 Mev , is assumed due primarily to the scarcity of data 

(particularly from in situ measurements) in Vette T s model. To 

compensate for this lack of data, the AE2 environment was averaged 

with the 0005 data to obtain a new model (referred to as the V-0 

model), which is more representative of the electron fluxes in the 

region 5 £ L £ 7 for geomagnet ically quiet periods (K < 1 ) and 

P 

solar cycle minimum conditions. The validity of this approach is 

supported by the close agreement between the differential energy 

spectrum obtained for L = 6 ,5 and the data given by Vette et al 

[1967] for L = 6.6. (Vette obtained results at this L shell from 

satellites in synchronous orbits.) 

The V-0 model is still not representative of the condition in 

the magnetosphere pertinent to this study, as the IMF-6 data were 

obtained predominantly during periods of sustained geomagnetic activity 

which resulted in 1 £ K £ 2. Furthermore, the radio noise data were 

P 

taken during March through October of 1971, a year which was charac- 
terized by an abnormally high amount of solar and geomagnetic activity 
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[see the ESSA publication. Solar- Geophysical Data, for that year]. 


Thus, we must extend the V-0 model on two accounts. 

a. First, it is necessary to compensate for differences 
in the average K p range represented by this model (K^ ^ 1 ) 
and the range from which IMP-6 data are obtained. The appro- 
priate factor by which the V-0 model must be augmented is derived 
from the results published by Pfitzer et al [1966], Freeman 
[1964], Owens and Frank [1968], and Craven [1966]. These authors 
show that there exists a positive correlation between an 

increase in K and the intensity of trapped electrons. In 
P 

particular, Craven and Freeman observed that, on the average, 
the increase in the integrated electron flux for E £ 40 kev, 
when K p goes from less than l’ to 2, is a factor of four to five. 

As was pointed out earlier, during the early recovery 
phase of a geomagnetic storm, the increase in the entire spec- 
trum of the trapped radiation is approximately constant and is 
given by the enhancement of the lower energy electrons . (This 
constant increase in the electron intensities, for all E, is 
also shown in Pfitzer's et al [1966] data.) Therefore, it 
shall be conservatively assumed that the average V-0 model 
must be increased by a factor of three to be. valid for 

1“ £ K 2. 

P 

b. Second, Vette [1971] points out that the expected 
enhancement of the maximum flux intensities at E = 500 kev is 
a factor of approximately four as the solar cycle (and hence 
solar activity) goes from minimum to maximum. The 500 kev 
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segment of the electron spectrum, however, appears to be the 
least sensitive to changes in geomagnetic activity [Frank et al, 
1964; Pfitzer et al 5 1966] . Thus, as stated previously, one 

is led to suspect that the increases in the lower, 40 kev £ E 
£ 200 kev, and higher, E £ 700 kev, electron intensities will 
be greater than four, as they respond much more dramatically 
to geomagnetic disturbances [Frank et al , 1964; Craven, 1966], 
Consequently, we shall again assume that the averaged V^O 
model must be increased by a second factor of three. 

The final model used to represent the quiet time (K < 1 ) 

P 

electron fluxes expected during 1971 is, therefore, the V-0 model 
increased by a factor of three. This model is then increased by a 
second factor of three in order to make it valid for fluxes at the 
geomagnetic equator during periods of geomagnetic activity resulting 

in 1 £ K £ 2. 

P 

3. Differential Electron Flux Spectrums for 3 £ K p £ 4 in the Geo^ 
magnetic Equatorial Plane 

The differential electron density spectrums for 3 £ K £ 4 are 

P 

derived from the results of Pfitzer et al [1966] and Frank et al 
[1964]. These studies show that the average increase in the electron 
intensities (for 40 kev < E < 2 Mev) from quiet time values to periods 
of moderate geomagnetic activity (K^ « 3) is 1.0 to 1.2 orders of 
magnitude. This same enhancement was also found by Craven [1966] 
and Freeman [1964] for the integrated electron fluxes, E ^ 40 kev, 
under similar changes in geomagnetic activity. 
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We have included in Figures 2.3a-b data obtained by 0G0-5 in 
the morning magnetosphere during a period of sustained geomagnetic 
disturbances in which K ^ « 3 . Comparing Figures 2.2b and 2.3b shows 
that the lower energy electron intensities are increased by a factor 
of ten; this is consistent with Pfitzer's and Frank 1 s studies. The 
higher energy electron population (E ^ 700 kev) , in turn, is enhanced 
by factors of decreasing magnitude, as both L and E increase. Com- 
bining these OGO-5 data with the published results , the average 
differential energy spectrums expected for 3 ^ K p £ 4 are taken to 
be equal to the model used in the preceding section, multiplied by 
a factor of eight. This model not only satisfies the sources quoted 
above, but also agrees with the average enhancements of trapped 
electrons caused by the storm time injection process summarized in 
the first paragraphs of Section II. B. 

An interesting feature of the data given in Figures 2.2a and 2.3a 
is the depletion of the fluxes near L = 3.0. The location of the 
minima in these figures is a function of the particle energy. 

However, for the ranges considered in this and the preceding sec- 
tion there is a consistent decrease in the electron fluxes (for all 
energies) in the region 2.5 £ L £ 4, compared to 1.5 £ L £ 2 and 
5.0 < L < 7. This segment of the magnetosphere, referred to as the 
slot region [Vette, 1971], serves to differentiate between the inner 
and outer radiation belts, and proves to be important to the discus- 
sion in the next section. 
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Figure 2.3 a-b DIFFERENTIAL ELECTRON FLUX SPECTRUMS FOR PARTICLES IN 
THE MORNING MAGNETOSPHERE DURING MODERATE GEOMAGNETIC 
ACTIVITY. MEASUREMENTS FOR 3<L<7 WERE MADE WHEN 
-8°<\ < 8° [From West, 1972]. 




4, Differential Electron Flux Spectrums for 5 £ K p £ 6 + in the Geomagnetic 
Equatorial Plane during the Early Recovery Phase of a Geomagnetic 
Storm 

Due to the scarcity of data for periods of sustained geomagnetic 

activity in which 5 £ £ 6 + during May - October 1971 (the time 

covered by the IMP-6 data), we shall develop a typical differential 

flux spectrum for the early recovery phase of geomagnetic storms which 

result in a K index of magnitude 5 £ K £ 6 + . This model will allow 
P P 

the comparison between theoretical calculations and the radio noise 
received by the IMP-6 during the 16-18 May and 24-28 September 1971 
geomagnetic disturbances. 

The storm time temporal variations of the electron intensities 
in the outer radiation belt follow a very intricate, but repetitive, 
pattern during the onset and main phase of large geomagnetic storms 
[Freeman, 1964; Frank et al, 1964; Craven, 1966; Owens and Frank, 

1968] , During the early recovery phase, the injection and depletion 
processes (which result in large, short term fluctuations) have 
generally subsided, and the fluxes reach their stable disturbed 
levels from which they gradually decay to pre-storm conditions 
[Freeman, 1964]. This early recovery phase is characterized by 
enhancements in the low, 40 kev < E < 200 kev , and high energy, 

E> 1.6 Mev electron fluxes by factors around 40 [Frank et al , 1964], 

Freeman [1964] saw this same general increase in intensities; his 
results for Kp ^ 4 show rises in these two components of the trapped 
electron population by factors of 50 to 100. Craven [1966] conducted 
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a detailed correlation between and the integrated electron density 
for E ^ 40 kev. His study showed that for Kp « 5 + , the increase was 
approximately 60 over quiet time values. 

Based on these results, we shall assume average electron fluxes 
which are a factor of 60 above the quiet time model developed for 
Kp < 1 . Even though this factor is probably not valid for electron 
fluxes of 200 kev ^ E £ 700 kev, the fact that an average model is 
used results in a calculated brightness for the magnetosphere which 
differs by less than a few percent from that obtained by using a 
model where factors of approximately 70 were taken for 40 kev < E 
£ 200 kev and E ^ 700 kev and a factor of 10 to 20 for 200 kev £ E ^ 
700 kev. 

One further important aspect of large geomagnetic storms on the 
trapped electron population is the injection of energetic electrons 
into the slot region. In this segment of the magnetosphere, increases 
in low energy electron fluxes by 1000 are common, and during the 
recovery phase similar increases (following a main phase depletion) 
are observed for the high energy electron intensities [West et al, 
1970]. This effect is clearly seen in Figures 2.4 and 2.5, which 
show 0G0-5 data taken during the geomagnetic storm of 11 June 1968, 
These data, which are typical of storms resulting in K p « 5 + [West, 
personal communication] , show that the injection process can bring 
flux levels in the slot region to values comparable with those in the 
outer radiation zone. This effect is only seen, however, for storms 
resulting in moderate to severe geomagnetic disturbances [Hess, 1968], 
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Figure 2.4 THE RESPONSE OF THE ELECTRON FLUXES IN THE SLOT REGION, 


L=3 . 5 , TO A MODERATE GEOMAGNETIC STORM. LABELS ON THE 
CURVES INDICATE DAYS PRIOR TO OR AFTER THE ONSET OF THE 
DISTURBANCE WHICH OCCURRED ON JUNE 11, 1968 [H. West, 
unpublished data, 1972]. 

Figure 2.5 THE RESPONSE OF THE ELECTRON FLUXES IN THE OUTER RADIATION 
BELT, L=4 . 5 , TO A MODERATE GEOMAGNETIC STORM. LABELS ON 
THE CURVES INDICATE DAYS PRIOR TO OR AFTER THE ONSET OF 
THE DISTURBANCE WHICH OCCURRED ON JUNE 11, 1968 [H. West, 
unpublished data, 1972]. 
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The complete model for the electron intensities existing during the 

early recovery phase of such geomagnetic disturbances must therefore 

be extended to include the fluxes to L ^ 3.1, These fluxes, for 

3,1 £ L £ 4.0 and 100 kev £ E £ 700 kev, were made a factor of three 

to four below the fluxes in the same energy interval at L & 5.0, and 

this factor was decreased monotonically , for increasing E, to unity 

for E> 1.0 Mev. This allowed the composite model for 5 £ K £ 6 + 

P 

and 3,1 ^ L £ 7.0 to be consistent with the data presented in Figures 
2.2b, 2.4 and 2.5 and to still retain the slot region as a separate 
entity from the outer radiation belt. 

5 . Pitch Angle Distributions 

The pitch angle distribution for the energetic particles proves 

to be an important consideration in the analysis to follow. Therefore, 

pitch angle data taken on the fluxes of Figure 2.2a are given in 

2.6a-c. These distributions are representative of particles drifting 

on closed L shells in the morning magnetosphere. It is possible, 

however, to have pitch angle distributions which are either isotropic 

or butterfly in nature [West, 1972], For this study, it will be 

shown qualitatively that the brightness of the magnetosphere is least 

for pitch angle distributions which are humped in form (i.e., the 

o 

maximum fluxes occurring at a = 90 and decreasing monotonically for 
O' < 90°), This distribution is a worst case condition and will be 
used in the quantitative portions of this paper. 

For L values 3.5 £ L £ 7, Figures 2.6a-c indicate that the 
particle fluxes at a given energy are within a factor of two to three 
for a « 45° compared to a = 90°. To put this observation into more 
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Figure 2.6 ELECTRON FLUX PITCH ANGLE DATA TAKEN CONCURRENTLY WITH 
THE DATA GIVEN IN FIGURES 2.3 a-b [From West, 1972]. 
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2 

quantitative terms, these data were approximated by a sin a law. 

That is, the flux of electrons at a given energy (E) and pitch angle 
(Qf) is found from the flux for O' = 90° by 

J (E) = J (E) sin 2 or (2.1) 

O' _L 

This approximation agrees very well with the data shown in these 
figures as well as the pitch angle data published by Katz [1966]. 

C . The Thermal Electron Population 

The effects of the thermal electron population on the calculated 
brightness of the magnetosphere have already been discussed in Chapter 
I, where it was shown that for regions in which the real part of the 
refractive indices, n , is less than or equal to zero, gyro-synchrotron 
radiation will not be generated. Furthermore, the radio model proved 
that if n + < 1.0, the intensity of the radiation at a given frequency 
is substantially reduced below the free-space situation. 

Once the electromagnetic energy has detached itself from the 
source, however, the effect of the thermal plasma reduces to the 
well-known problem of electromagnetic waves propagating through an 
ionized medium [Budden, 1961; Ginzburg, 1961; Ratcliffe, 1962; Stix, 
1962]. The portion of this theory that is relevant to this study 
pertains to the propagation of radio waves at a frequency through 
regions where n^f^) 0. It is shown by these authors that if the 

conditions shown in Figure 2.7 exist, the observer will not receive 
the energy incident on the left side of the region (assuming the width 
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REGION I 
n ± (f,) >o 


REGION 2 


REGION 3 


n±(f,)£o 


PROPAGATING ELECTROMAGNETIC 
WAVE OF FREQUENCY f, 


OBSERVER 


DIRECTION OF DECREASING n+ 


W 


Figure 2.7 RADIATION OF FREQUENCY LESS THAN OR EQUAL TO f IN REGION 
ONE IS SHIELDED FROM THE OBSERVER BY REGION TWO WHERE 
n + (f < f ) = 0. 

(w) is many wavelengths of the incident radiation). For an observer 
located in the interplanetary medium, the situation depicted in this 
figure will occur for radio noise, 30 kHz £ f ^ 300 kHz, generated 
within certain segments of the magnetosphere (as discussed in Section 
II. G) . Thus, it is necessary to determine the thermal electron 
density in (B,L) space as a function of geomagnetic activity, as this 
particle population will be significant to the results of the next 
chapter. To this end, in situ and indirect measurements on the 
density and response of the ambient plasma to magnetic disturbances 
are used to derive an appropriate model. 
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1 . Response of the Thermal Plasma at the Geomagnetic Equator to 
Variations in Kp 

The plasmapause, defined in the literature as the location of a 
sudden decrease in the thermal plasma density as one moves radially 
outward from the earth, separates a region of high thermal electron 
density (the plasmasphere) from a region of very low density (the 
plasraatrough) . (Schematic diagrams defining these regions in the 
earth's magnetosphere are given in Chappell [1972] and Carpenter [1968].) 
This phenomenon has been studied to a great extent in the geomagnetic 
equatorial plane, and the results are briefly summarized below: 

a. 0G0-1 and 0G0-3 data [Taylor et al . , 1968], obtained for 

local time near midnight, show a sudden decrease in the magneto- 

9 “3 6 “3 

spheric proton density from 10 ions m to 5 x 10 ions m in 

instances of less than one half an earth radius as one moves 

radially outward from the earth along the geomagnetic equator. 

6 -3 

(The lower bound of 5 x 10 ions m is a result of instrument 
limitations.) Furthermore, there is a correlation between geo- 
magnetic activity and the distance to the plasmapause, which 
decreases from L = 5.0 during undisturbed times to L = 3,5 or 
less for Kp » 5 + . 

b. Whistler data [Carpenter, 1966, 1970, 1971; Angerami 
and Carpenter, 1966; Carpenter et al . , 1971] show these same 
results with the added feature that a bulge in the plasmapause 
exists for local times between 1800 and 2100 hours. At this 
bulge, the distance to the plasmapause may be as large as 5.0 
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earth radii, compared to 3.5 for the 0600 - 0800 time zone, 
during periods of moderate geomagnetic activity. Whistler data 
also indicated that the position of the plasmapause changes with 
a delay of at most a few hours in response to geomagnetic dis- 
turbances . 

c. A comprehensive study of the plasmapause has been 
carried out recently by Chappell et al [1970] and Chappell 
[1972] from OGO-5 data. They showed that the response of the 
plasmapause to changes in K p is observed both in the 1000 ± 2 
hour and 0200 t 2 hour local time zones after a two to six hour 
time delay from the onset of the geomagnetic disturbance. This 
time delay was most evident for the 0200 + 2 hour zone, but 
alsofits the data from the 1000+ 2 hour region. 

As a result of the increased sensitivity of the experiments 
aboard the OGO-5 (compared to OGO-1 and 0G0-3) it was also 
discovered that , for undisturbed periods , the change in the 

thermal proton density (as the plasmapause was crossed) was 

9 —3 5 —3 

from 10 ions m to approximately 3 x 10 ions m over a 

distance of one earth radius. As geomagnetic activity increased, 

this same three orders of magnitude change in density occurred 

over distances of only 600 km. 

For this paper, the results of these investigations were 
combined to obtain the following relationship between the location 
of the plasmapause in the geomagnetic equatorial plane and K^. 

Lp — 5.6- 0.45 K p (2.2) 
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We shall assume (based on OGO-5 data) that, after a maximum delay of 
six hours, Equation (2,2) gives the average plasmapause location in 
the 2100 to 1800 hours local zone. From 1800 to 2100 hours, it will 
be assumed there exists a bulge in the plasmapause extending to 
L « 5.0. 

The data obtained by OGO-5 were also used to determine density 
profiles for the thermal electron population as a function of Kp. 
Linear approximations to these data, for the three Kp ranges specified 
previously, are given in Figure 2.8 

D. The Electron Plasma and Gyrof requencies as a Function of Kp and L 
in the Geomagnetic Equatorial Plane 
The electron plasma frequency, as a function of Kp, is found from 
the data in Figure 2.8 by using [Rossi and Olbert , 1970] 

f « 9 / N Hz (2.3) 

P N e 

where 

-3 

N = electrons m 
e 

The results are shown in Figures 2.9 a-c. 

To determine the electron gyrof requency at the geomagnetic 
equator, the variation of the magnetic field intensity (B) in this 
plane as a function of K p must be obtained. Roederer [1970] states 
that the field at the magnetic equator can be expressed as 

k o 

B 0 (r) = — + k^ ~ k 9 r cos 0 (2.4) 
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Figure 2.8 THERMAL ELECTRON DENSITY IN THE GEOMAGNETIC EQUATORIAL 

PLANE AS A FUNCTION OF L FOR THREE K RANGES [From Chappell 

P 

et al , 1970], 
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where 


4 

= 3,11 x 10 gamma 

12 x 10 3 < * 1 ‘° gamma = 10 gauss) 

k gamma 

1 R J 

S 

, 2.27 x 10 4 

k 2 = gamma 

R s 

R s ~ distance from the center of the earth to the subsolar point 
of the magnetosphere in earth radii. 

0 = longitude east of midnight. 

r = distance from the center of the earth at which B is to be 

o 

evaluated (in earth radii). 

Changes in Kp are accounted for in this equation by R s . For K p £ 2 , 

R g « 10 R e ; during a severe compression of the magnetosphere, R may 

s 

be as low as 6 R fi (earth radii) . 

This equation was evaluated for R c = 10 R n , K <, 2: R = 8 R„ 

JJi p 9 s E 9 

3 ^ Kp ^ 4; and R s = 7 Rg, 5 £ £ 6 + as a function of 0 . The 

values for B Q were then used to calculate the electron gyrof requencies 
for nonrelat ivistic particles in the equatorial plane; the results 
are plotted in Figures 2.9a-c with the electron plasma frequencies. 

The curves for and f allow an important simplification to 
be made. Intuitively, we expect the radiation emitted by super thermal 
electrons to increase for decreasing L in the geomagnetic equatorial 
plane, since both the magnetic flux density (B) and the energetic 
electron population increase. This indeed occurs, as demonstrated in 
Chapter III, where the calculations show that the bulk of the radiation, 
for any K p value, emanates within two to three earth radii from the 
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Figure 2.9 a-c ELECTRON PLASMA AND GYROFREQUENC IES IN THE GEOMAGNETIC EQUATORIAL PLANE AS A 

FUNCTION OF L FOR THREE K RANGES. THE UPPER HYBRID FREQUENCY (f ) WAS CALCULATED 
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USING f AND f (DIPOLE). 
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plasmapause (in a direction away from the earth). The preceding 
figures illustrate that, for L p + 3 R the difference between the 
electron gyrof requency as one proceeds from the midnight to the noon 
meridian is a maximum of 3.5 kHz (for each Kp range) with respect 
to the dipole field. Since this discrepancy proves insignificant to 
the theoretical results, the earth's magnetic field can be satisfac- 
torily approximated by a centered dipole model for this study. 

E. Extension of the Particle Models to Regions off the Geomagnetic 
Equator 

Thus far we have developed the models for the thermal and 

energetic electron populations as a function of K in the geomagnetic 

P 

plane. These data, however, are not sufficient to calculate the 
brightness of the magnetosphere. To complete these models, the dis- 
cussion of the previous section must be extended to regions of the 
magnetosphere off the equator. 

1 . The Thermal Electron Density in (B.L) Space 

Carpenter [1966] has found that the plasmasphere can be modeled 

as a doughnut-like cloud of ionization with a sharp, field-aligned 

boundary. This magnetic field line, which bounds the region of high 

thermal electron density, passes through the plasmapause at L in the 

P 

geomagnetic equator. 

This boundary serves as a natural divider for calculating the 
thermal electron densities off the geomagnetic equator. Well within 
the plasmasphere, the thermal electron density is high enough to cause 
collisions between electrons to establish a diffusive- 
equilibrium distribution of plasma along the field lines. Outside 
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the plasmapause, collisions are rare, and an idealized collisionless 
distribution is appropriate. These features imply that, for magnetic 
field lines within the plasmapause and near the earth's surf ace, the 
plasma density decreases with distance along the field line (measured 
from the earth) in accordance with a diffusive distribution law; 

-4 

outside the plasmasphere , the density decreases with distance as R 
[Park, 1970; Helliwell, 1970], As the portion of the magnetosphere 
outside of the plasmapause proves important to this study, the thermal 

electron density in (B,L) will be calculated by extrapolating the 

-4 

value at the equator (for a given L) using this R power law. 

2 . The Energetic Electron Population in (B,L) Space 

In this paper, only energetic particles which are stably trapped 
will be considered sources of gyro- synchrotron noise. The lifetime 
of these particles, in the stable trapping regions, ranges from three 
days at L « 6 to six days at L « 4 [Hess, 1968]. Consequently, the 
energetic electron population would be lost if it were not for the 
continual replenishment of these particles by geomagnetic disturbances. 

As a result of the short bounce and drift periods for electrons 
with 100 kev £ E £ 5 Mev (which are a maximum of seconds to hundreds 
of minutes, respect ively, [Hess , 1968]) and the long particle lifetimes, 
these particles will distribute themselves uniformly along drift shells 
within hours after injection. The resulting population in (B,L) 
space can therefore be determined by using Liouville's theorem, which 
is given by 
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J p (E t COS CV p ) = Jq (E , cos = jJe, 


1 f- (1 - cos 2 
B P 




(2.5) 


This equation shows that the flux of particles, J p (E, cos a), having 
energy E and pitch angle or measured at a point P on a given magnetic 
field line, results from the flux of particles, J Q (E, cos ) on the 
same field line, with pitch angle cy f passing through the point O. 
Equation (2.5) can be put into the following more useful form when 
a p = 90° 


Jp(E.O) 




( 2 . 6 ) 


where, as a result of the dipole field approximation , 


b p = 


B q (1 + 3 sin 2 X m )2 


6 

cos X 


m 


(2.7) 


In these equations 

Bq = magnetic flux density at the geomagnetic equator measured 
at point 0. 

Bp = magnetic flux density measured at point P along the same 

dipole field line passing through point O. 

X = geomagnetic latitude 
m 

Jp(E,0) = perpendicular flux of particles of energy E measured 
at the point P. 
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B„1i] 


1 - 


= flux of particles of energy E measured on the 
equator (at point 0) with pitch angle a defined 


by cos a = II 


b a 




For example, if we consider L = 5.0 and £ 2 , Figure 2.2a 

, 7 -2 -1 -1 -1 

shows that Jj_(E = 266 kev) « 4 x 10 electrons m kev ster sec 


at \ « -15 , To calculate the perpendicular flux at L = 5.0 and 


m 


1=3 8 we proceed as follows: 

m * 


Jp( E , cos a p ) = J P (E,0) = J p (E,0,L = 5°, = 38°) 

- V E - cos V L = 5 > K - °°> - J 0 | K - 

= J 0 [e,(1 - 0.16)*] = J 0 (E,0.91) 

2 

Using the sin a pitch angle distribution, 

J o (E,0. 91) = J 0 (E,0) sin 2 23.9° = 0.16 J Q (E,0) 

Going through this procedure once more to account for the flux of 

Figure 2.2a being measured at A = -15°, we find that: 

m 

J(E,0,L = 5.0, X = 38°) = (0.25) J(E,0,L = 5.0, \ = -15°) 

m m 

,7 -2 - 1,-1 -1 

= 10 electrons m sec kev ster 

This value compares well with the 0G0-5 data shown in Figure 2.10, 
which is the extension of Figure 2.2a to higher geomagnetic latitudes. 
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Figure 2.10 EXTENSION OF THE ELECTRON FLUXES OF FIGURE 2.2a TO 

REGIONS OFF THE GEOMAGNETIC EQUATORIAL PLANE [H. West, 
unpublished data, 1972]. 




F. Stratification of the Magnetosphere 

To evaluate the equation for the transfer of radiation (1,11) 
from source to observer on a digital computer, it is necessary to 
stratify the magnetosphere. Previous sections pointed out that a 
sufficiently accurate value for the generated radio noise can be 
obtained by considering only the electrons trapped in a dipole magnetic 
field. The particle dynamics in this type of field configuration have 
been investigated by Roederer [1967, 1970], who showed that the 
trapped particles are constrained to move along closed L shells (or 
drift shells) which are determined by rotating dipole field lines 
about the dipole axis. Therefore, the most natural stratified repre- 
sentation of the magnetosphere is a collection of concentric L shells. 
It remains, however, to determine a convenient distance (Al) between 
each shell at the geomagnetic equator. 

This problem was resolved by running a few trial calculations 
for the noise based on different AL values, under the assumption that 
all the radiation emanates from the geomagnetic equatorial plane. The 
results indicated that if AL was greater than 0.3 earth radii, the 
calculated radio noise spectrum varied considerably from one shell to 
the next, and the summed spectrum for all the shells showed a charac- 
teristic cyclotron harmonic structure at the lower frequencies. This 

is a result of these frequencies being close to the fundamental 
gyrof requency of electrons with E « 100 kev. [See Bekefi, 1966; 

Figure 6.10a, for an example of this effect.] The cyclotron struc- 
ture was particularly evident for 5 £ ^ 6 + , as f 0 attains values 

to 66% of the lowest frequency considered in the paper (30 kHz). 
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However, these fluctuations in the summed radio spectrums become 
insignificant and a reasonably good power law for the radio noise was 
obtained when a A L value of 0.125 earth radii was used. The magneto- 
sphere was therefore modeled as a set of concentric drift shells 
separated by this AL value in the geomagnetic equatorial plane. The 
center of each shell in this plane (in the region of interest) is 
given by 

L s = 3.0 + j (0 . 125) j - 0,1,2 ... (2.8) 

G. The Refractive Indices in (B,L) Space 

Using Equation (1,38), the values for n + (both when theta equals 
zero and theta equals ninety degrees) were calculated and plotted in 
Figures 2.11 a-f. These curves show that, for each mode, there exists a 
region at the geomagnetic equator for which n + is equal to zero. Further- 
more, for each K p range, this region is located within the plasmasphere 
and ends abruptly at the plasmapause. This same situation also exists 
when the refractive indices are calculated in (B,L) space along the 
plasmasphere- plasmat rough boundary, Heurist ically , this result is sup- 
ported by the following argument. Choose the L value (for a given K 

P 

range) which determines the location of the plasmapause in the geomag- 
netic equatorial plane. Moving toward the magnetic poles along field 
lines which pass through this L value, both the field flux density (Equa- 
tion 2.4) and the thermal electron population increase. Thus, at a given 
frequency, the real parts (n + ) of the refractive indices will decrease in 
magnitude toward zero for a value less than or equal to unity at the equa- 
tor. If is already zero at the equator, it will not change as one moves 
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REFRACTIVE INDEX (n + ) 


L>4.75 


L >4.75 




(a) (b) 

Figure 2.11 a-f REFRACTIVE INDICES (n ) VERSUS FREQUENCY, L, K AND THETA IN THE GEOMAGNETIC 

+ p 


EQUATORIAL PLANE. 







REFRACTIVE INDEX (n + ) 




(c) (d) 


Figure 2.11 (continued) 









along the field line. (This argument is based on the assumption that 
the frequency of observation is always greater than the upper hybrid 
frequency encountered along the field lines defining the plasmapause. 
This assumption will always be valid for regions in the magnetosphere 
from which appreciable radiation will be shown to emanate.) Conse- 
quently, for frequencies between 30 kHz and 300 kHz, radiation from 
within the plasmasphere will not reach an observer located in the 
interplanetary medium. The noise that would be received, therefore, 
comes from energetic particles located beyond the magnetic field line 
defining the plasmapause (i.e., primarily from the outer radiation belt 
electrons since the inner belt, located between L = 1.1 and L « 2.0, 
is found within the plasmasphere for K p < 6 + : see Equation 2.2). 

Figures 2.9 a-c, 2.10 a-c, and 2,11 a-c also show that the 

refractive indices (for 30 kHz ^ f ^ 300 kHz) vary monotonically from 

zero to one over a very short distance as the plasmapause is traversed 

from the plasmasphere to the plasmatrough . The models of Chapter I 

proved that, for regions where n < 1.0, the power radiated by an 

electron is decreased below that which would be emitted if n = 1 0 

± 

It is therefore possible to limit this study to regions where > o.l 
without introducing significant error in the final result. This state- 
ment is supported as follows. The segment of the magnetosphere, in 
which 0 S n± £ 0.1, has a width of at most 0.125 earth radii (measured 
in the geomagnetic equatorial plane) while the overall emission region 
(where by definition 0.1 ^ n ^ 1.0) has an extent of at least 2.5 
earth radii. Thus, the calculations for the next chapter (see Section 
III.B) will show that the noise generated within the segment of the 
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magnetosphere adjacent to the plasmapause, is almost three orders of 
magnitude below that generated in the overall emission region; it can 
therefore be safely ignored. 

H. Ray Path Calculations 

For the frequency range 30 kHz to 300 kHz, the magnetospheric 
models valid for this study, and radio noise propagating in a region 
characterized by 0.1 £ n ± ^ 1.0, paragraph eight of I.B states that 
the transfer of radiation from the source to an observer can be 
calculated using the principles of geometric optics. That is, the 
magnetospheric noise received by an observer in the interplanetary 
medium can be calculated by determining the brightness received along 
ray paths that pass through the emission region and terminate at the 
observation point (see Figure 1.3). The calculations for these ray 
paths (referred to as raytracing) are greatly simplified by the 
following. 

1. The ray path geometry is easier to compute if we assume the 
paths start at the observer (i.e., the observer is modeled as a trans- 
mitter) and pass through the magnetosphere. By the reciprocity 
theorem, these same paths would be obtained if they started in the 
magnetosphere and were traced to the observer. This latter approach, 
however, complicates the calculations unnecessarily. 

2. Consider the ray paths of Figure 2.12, which are drawn in the 
geomagnetic equatorial plane. For points along the second path, before 
it hits the plasmapause, the refractive indices are approximately equal 
to unity (as shown in the preceding section). Therefore, the ray paths 
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Figure 2.12 RAY PATH GEOMETRY IN THE GEOMAGNETIC EQUATORIAL PLANE FOR 
AN OBSERVER WHOSE GEOCENTRIC DISTANCE IS 32 EARTH RADII 
(not drawn to scale) . 



Figure 2.13 EXPANDED VERSION OF FIGURE 2.12. THIS FIGURE DEFINES 
THE RAY PATH GEOMETRY PERTINENT TO THIS STUDY. 
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can be taken as a straight line from the observer to approximately one 
L shell from the plasmapause. 

3. As the ray approaches the plasmasphere , it encounters a medium 
in which the electron density increases rapidly in a direction toward 
the center of the earth. Furthermore, the rate of change of the 
electron density (and hence plasma frequency) is much greater than the 
rate of change of the electron gyrof requency . This results in a 
decrease in magnitude for n , which is primarily due to the electron 
density gradient. Consequently, the magnetic field can be neglected 

in the ray path calculations. 

4. Consider the geometry shown in Figure 2.13. Define the L 
shell which contains the plasmapause as (L shell) „ For a ray path 
hitting the plasmapause, the distance it penetrates the plasmasphere 
is designated as , and the distance from the ray entry point to the 

exit point from (L shell) is X . Provided both X and Z are small 

m mm 

compared to the distance (D^) defined in the figure, the magnetosphere 
can be modeled as a plane (at the ray entry point) in which the ther- 
mal electron density increases linearly, with decreasing L, along the 
earth-observer line. Based on this model, it is possible to obtain 
explicit solutions for X m and Z^ [Kelso, 1964] which are given by 

f 2 cos 2 ¥ (2.9) 

Z m TU2 

4f 2 sin 2 ¥ cos ¥ (2.10) 

x 

m T| A 2 
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where 


f = frequency of the electromagnetic wave propagating along the 

ray path (Hz) 

2 3-9 

A “80.6 meters sec ~ 

T| = slope of the electron density profile (electrons m 
¥ = angle of incidence of the ray path as defined in Figure 
2.13. 

Using these equations, and were calculated for ray paths starting 

at the observer and passing through the three model magnetospheres 

developed in this chapter. In the worst case, where 2^ assumed its 

largest value, it was found that Z m « 640 km. The value for X , how- 
9 m m 9 

ever, increases as the ray path becomes tangent to (L shell) 

(i.e., ¥-*90°). If this case is neglected, by setting ¥ ^ 80°, the 
worst case results in 3^ » 2800 km. Since the corresponding value for 
D m is approximately 32,000 km, the approximations used in calculating 
the ray paths (for ¥ ^ 80°) are justified. 

When 80° < ¥ < 90°, trial calculations showed that the values 
obtained for 3^ by Equation (2.10) were at most 10% less than those 
which would actually be obtained if a rigorous raytracing procedure 
were used. If this 10% discrepancy is accounted for in the next 
chapter, however, a change in calculated brightness of the magneto- 
sphere of less than 1.0% would be obtained. This reduction in the per- 
cent error is due to the fact that 3^ is always much less than the 
overall length of the ray path in the emission region; hence the 
integration of the emission coefficient J(i) along each ray path 
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reduces the relative magnitude of this error. Since the 1.0% error in 
the brightness calculations is insignificant, the approximations for 
Z m and given by Equations (2.9) and (2,10) were used for obtaining 
all the ray paths. 

Examples of the calculated ray paths (drawn to scale) are illus- 
trated in Figures 2.14, 2.15 and 2,16 for the three K p ranges discussed 
previously and for a frequency of 110 kHz. In each figure, two sets of 
drawings are given. One represents ray paths in the geomagnetic 
equatorial plane, and the second represents ray paths in the plane 
defined by the earth's dipole axis and a line drawn from the center of 
the earth to the observer. It should be noted that in this second set 
of drawings, a number of ray paths (particularly for 5 £ K p £ 6 + ) were 
terminated before they traversed the entire magnetosphere. Such paths 
are those which come close to the earth's surface at large geomagnetic 
latitudes. They were truncated because they start to penetrate regions 
of high thermal electron densities and increased magnetic field 
strengths, which invalidate the simplifications listed above. Their 
termination, however, is not of major concern since, if these ray 
paths are extended, they continue through portions of the magnetosphere 
where the trapped energetic electron population is orders of magnitude 
less than that which exists along the segments of the paths drawn in 
these figures. The truncation of these paths will therefore introduce 
second order errors in the final results obtained in the next chapter. 
In fact, calculations in the next chapter will show that all ray paths 
could be terminated at a geomagnetic latitude of ~ ± 40° (resulting 
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Figure 2.14 RAY PATHS CALCULATED IN THE 0-A PLANE (DEFINED IN THE 
TEXT) AND THE GEOMAGNETIC EQUATORIAL PLANE FOR GEOMAG- 
NETICALLY QUIET DAYS (drawn to scale), cr and € are the 
observer look angle, subtended between the launch direc- 
tion of a particular ray and the earth-observer line, 
projected onto the 0-A plane and the geomagnetic plane 
respectively . 
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Figure 2*15 RAY PATHS CALCULATED IN THE 0-A PLANE (DEFINED IN THE 
TEXT) AND THE GEOMAGNETIC EQUATORIAL PLANE DURING SUS- 
TAINED GEOMAGNETIC ACTIVITY RESULTING IN 3 < K <4 

P 

(drawn to scale), or and e are the observer look angle, 
subtended between the launch direction of a particular 
ray and the earth-observer line, projected onto the 
O-A plane and the geomagnetic plane respectively. 


97 



Figure 2.16 RAY PATHS CALCULATED IN THE 0-A PLANE (DEFINED IN THE TEXT) 

AND THE GEOMAGNETIC EQUATORIAL PLANE DURING GEOMAGNETIC 

ACTIVITY RESULTING IN 5 < K < 6 + (drawn to scale) . 

P 

or and € are the observer look angle, subtended between 
the launch direction of a particular ray and the earth- 
observer line, projected onto the O-A plane and the 
geomagnetic plane respectively. 
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CHAPTER III 


THEORETICAL CALCULATIONS 


A. Introduction 

Using the models given in the last two chapters, it is now 
possible to calculate the expected brightness of the earth’s outer 
radiation belt due to gyro- synchrotron radiation. These calculations 
are made only for the extraordinary mode, since previous results showed 
that the power radiated into this mode is, in general, greater than or 
equal to that radiated into the ordinary mode, which will therefore 
not be discussed in this chapter. 

Sections III.B through III.F neglect the bulge in the plasmapause 
and consider the case of an observer and ray paths located in the geo- 
magnetic equatorial plane. The complete problem, in which these two 
restrictions are relaxed, is discussed in Section III.G, and the final 
results are then presented in III.H. 

B. The Emission Coefficient, J(f) 

The emission coefficient, as a function of Kp, was calculated 
using Equation (1.57) and the results are plotted in Figures 3,1 a-c. 
These figures show that, as the plasmapause is approached, radiation 
at the lower frequencies is greatly reduced from that which exists in 
the plasmatrough. This decrease in the emission is a result of the 
Razin effect discussed in Chapter I (i.e., n ± (f)— ^0, for all frequen- 
cies, as L Lp) . It is interesting that the distance to the plasma- 
pause, from the L shell in which J(f) begins to decrease, is directly 
related to the magnitude of K . This phenomenon is a result of the 
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EMISSION COEFFICIENT (woMs/m 3 -Hz-sier) 



Figure 3.1 a-c THE EMISSION COEFFICIENT J(f) IN THE GEOMAGNETIC 


EQUATORIAL PLANE AS A FUNCTION OF FREQUENCY FOR 
VARIOUS L AND K VALUES. 





increase in the thermal electron density gradient as increases 
(see Figure 2.8). That is, as the density gradient increases, n ± (f) 
goes from unity to zero over a shorter distance. The maximum value for 
this distance (measured from the plasmapause to the L shell in which 
the Razin effect, at any frequency, is negligible) occurs for 
1~ :£ Kp £ 2, and is equal to approximately AL (AL = 0.125 earth radii). 

C. Path Lengths as a Function of Look Angle, e 

Figure 3.2 gives the length of the ray paths, within each L shell, 

as a function of the observer look angle (e). This angle, shown in 

Figures 1.5 and 2.18, is defined as the one between the observation 

direction and the earth-observer line (i.e p , each value of e defines 

a particular ray path in the geomagnetic equatorial plane). The sudden 

increase in the path length, in each L shell, occurs when a ray passes 

tangential to that shell. This takes place, for an L shell given by 
* 

L = L , when e has the value 

L* - 0.0625 ^ (3.1> 

* J 

where 



D = distance from the center of the earth to the observer 
o 

located in the geomagnetic equatorial plane, measured in 
earth radii. (D 0 « 32 earth radii in this study.) 
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LOOK ANGLE « (degrees) 

Figure 3.2 PATH LENGTHS IN THE GEOMAGNETIC EQUATORIAL PLANE (WITHIN A 
GIVEN L SHELL) AS A FUNCTION OF THE OBSERVER LOOK ANGLE. 


As an example, for K 3 , the ray length in the L = 4.0 shell is 0.25 

P 

earth radii when e = 0°, but increases to 2.0 earth radii for e « 7° . 
When Kp < 3, however, the path length in this shell is zero since it is 
located within the plasmasphere, and ray paths for 30 kHz £ f £ 300 kHz 
were shown not to penetrate this portion of the magnetosphere. Thus, 
only path lengths in L shells, where L is greater than L^ (given by 
Equation (2.2)), have values other than zero. 
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This fact is shown in Figure 3.2 by the arrows, which indicate 

the location of the plasmapause as a function of K . For a given K 

P ^ 

range, shells to the left of each arrow are shielded from the observer, 
and hence do not contribute to the total length of any ray path that is 
calcul ated, 

D. The Absorption Coefficient, K(f) 

The next quantity to be considered is the absorption coefficient, 
which was calculated from Equation (1.61) for f ~ 30 kHz and the three 
Kp ranges previously defined. Since K(f) is proportional to the 
emission coefficient, which increases when the frequency decreases 
(see Figures 3.1 a-c) , and is inversely proportional to the frequency 
squared, the choice of f = 30 kHz results in the maximum expected 
absorption of the radio noise for 30 kHz £ f ^ 300 kHz. Therefore, 

K(f ) = 30 kHz can be considered a worst case situation. Figure 3.3 
gives the absorption, in dB per earth radius, for this frequency. 

The dip in these absorption curves at L « 4.5 is a result of the 
slot region, which is located in the magnetosphere from 2.5 < L < 4.5. 
Within this region, the decrease in the low energy electron fluxes 
(compared to L ^ 5.0) causes a decrease in the emission coefficient 
j(f) at low frequencies. (This result will be discussed in detail 
in Section III.F.) Due to the proportionality between K(f) and J(f), 
a similar reduction in K(f) is also observed. The dip in K(f) at 
L w 4.5, however, should not be confused with the abrupt cutoff in 
J(f), and hence K(f), caused by the Razin effect at low frequencies. 
(The subsequent rise in the absorption, when L — > L^, is due to the 
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Figure 3.3 THE ABSORPTION COEFFICIENT K(f) IN THE GEOMAGNETIC 

EQUATORIAL PLANE AS A FUNCTION OF L FOR f=30 kHz AND 

VARIOUS K VALUES, 

P 


increase in the magnetic flux density (B) , which causes a corresponding 
increase in J(f) and K(f)). 

The results for K(f) in Figure 3.3 allow an important simplifica- 
tion to be made. Since the longest path length within any L shell is 
not more than 4 earth radii, the maximum possible attenuation suffered 
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by the radio noise propagating along any ray path (for f ^ 30 kHz) is 
- 8 

pa 4.0 x 10 dB. The gyro- synchrotron self- absorpt ion process, there- 
fore, results in a decrease in the radiation of such small magnitude 
that it can be ignored in this study. 


E , Brightness of the Outer Radiation Belt as a Function of 

The brightness observed along ray paths passing through the 
magnetosphere is derived from Equation (1.16). This equation, however, 
can be simplified on two accounts. First, the results of the last 
section showed that the absorption coefficient could be neglected 
(i.e., the optical depth (t) is equal to zero). Second, the emission 
coefficient is appreciable only for L shells given by L > L + 0.125, 

r 

2 2 

where calculations give 0.1 ^ n 1.0 and n » n, . Thus Equation 


(1.16) reduces to 


1(f) 



J(M) 
n ±(f M 


<U = 


I 


J(f ,1) 

-±< f ■» 


AL 


(3.2) 


This expression was evaluated as a function of look angle (e); the 
results for f b 30 kHz, 92 kHz and 300 kHz are given in Figures 
3*4 a-c. These figures illustrate the following general characteris- 
tics for 1(f), 

1. A peak in the calculated brightness occurs for the path 
passing tangential to the L shell given by + 0.0625 (i.e., this 

path is defined by = sin [L^/D^]), and the maximum is at least 
6 dB greater than the brightness observed along the paths defined 
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BRIGHTNESS (watls/nZ-Hz-ster) 



(a) (*>) < c > 

Figure 3.4 a-c BRIGHTNESS OBSERVED ALONG RAY PATHS PASSING THROUGH THE EARTH’S MAGNETOSPHERE, IN THE 
GEOMAGNETIC EQUATORIAL PLANE. EACH VALUE OF 6 SPECIFIES A PARTICULAR RAY PATH IN THIS 


PLANE. 






by e = 0° . This maximum in 1(f) is a result of the following. For 
e < e 1# the ray lengths within each L shell are less than the corre- 
sponding lengths for the path given by e . For paths given by e > , 

the number of shells traversed by these rays is less than that traversed 
by the path given by e = . Thus, the path defined by e = has the 

longest lengths within each shell while still traversing the maximum 
number of L shells, and the maximum brightness is consequently observed 
in its direction. 

2. For the lower frequencies and K £3, there exists a dip in the 

P 

calculated brightness at L « L +0,25. This decrease in 1(f) can be 

P 

explained by considering the trapped electron model of Section II. B. 3. 

In this model, the particle fluxes for 100 kev <. E ^ 700 kev within the 
slot region, are three to four times less than the corresponding fluxes 
at L « 5.0, while the higher energy fluxes are the same in both 
regions. Now, as shown in the next section, low frequency noise is 
generated predominantly by low energy electrons, and the high frequency 
noise (f > 92 kHz) is generated by high energy electrons (E > 600 kev). 
Thus, consider two ray paths, one passing through the slot region 
(L ps 4.5), the second passing through the heart of the outer electron 
belt (L w 5.0), From the above discussion it is apparent that the 
second path will have a greater brightness at low frequencies than the 
first, but both will be of comparable brightness at the higher 
frequencies (i.e., the dip does not appear for these frequencies). 

The decrease and subsequent increase (as L increases) in the bright- 
ness curves for f = 30 kHz is a consequence of the electron fluxes in 
the slot region, which emit detectable radiation when K >3, (Note 
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that this effect is not observed for 1 ^ K ^ 2 since for this case 

P 

the slot region is located within the plasmasphere . ) The relative size 

of the dip is also positively correlated to K^, which is due to a larger 

portion of the slot region being exposed as K increases. 

P 

3. Define the half power angular extent of the emission region, 

at a frequency f, by e (f), where I(e ,f) = 0.5[l(e = 0°,f)]. 

max max 

Assume also that for e > e , I(e,f) =0.0. (The definition and the 

max 

assumption simplify the analysis to follow and introduce no more than 
10% error in the final results.) Then, Figures 3.4 a-c show that 
e mQV (f) decreases slightly as f increases. That is, the emission 
regions are confined to segments of the magnetosphere whose outer 
boundary is progressively closer to the earth as the observation 
frequency increases . 

This relationship between e (f) and f is due to a monotonic 
increase in [ d[log I(f)]/d[log f ] | as both f and L increase. The 
increase in |d[log I(f)]/d[log f]|, in turn, is a consequence of the 
decrease in the magnetic flux density and the softening of the electron 
flux spectrums as one moves away from the earth (see Figure 2.2b). 

F. The Range of Electron Energies and Harmonic Numbers which Contribute 
Radiation at a Given Observation Frequency 

Equation (1.56) shows that, within a given L shell, the emission 
coefficient J(f) is obtained by performing a summation over harmonic 
numbers (s) . That is, to calculate J(f), s is stepped from unity to 
’infinity’, and the radiation from electrons, whose energies are 
determined by the harmonic number and frequency of observation, is 
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summed. (Recall the discussion in Section I.E on the Dirac delta func- 
tion in this equation. This function implicitly gives the energy of the 
electrons which radiate at a frequency f once s, the positional geometry 
of the observer relative to the earth, and an L shell are specified.) 

It is therefore possible to calculate the ratio, for each L shell, of 
the power radiated at a frequency (f) and harmonic (s) (and hence from 
electrons of discrete energy (E)) to the total emission coefficient 
J(f) for that shell. The results of these calculations, expressed as 
a percentage of J(f), are given in Figures 3.5 a-c and 3.6 a-c for 

3 <; K <; 4. 

p 

From the discussion of the preceding paragraph, it is evident that 
the plots given in these figures should not be continuous, but should 
have data only at integral values of s and discrete values of E, The 
intention of this section, however, is to show clearly the relationship 
between J(f) and these two quantities. Therefore, the envelopes of 
these discrete plots are given. 

These figures illustrate that within each L shell the emission 
coefficient at low frequencies results from low harmonic numbers and 
electron energies. Furthermore, as the frequency of observation (f) 
increases, both the range of harmonic numbers and electron energies 
which account for the radiation at f also increase. 

Thus, consider the ray path defined by e = 0° in the geomagnetic 
equatorial plane. The brightness 1(f) observed along this path 
is obtained by summing J(f) for each L shell the ray traverses and 
multiplying the total by a constant. From the results given in 
Figures 3.5a and 3,6a it can be concluded that for f = 30 kHz, 70% 
of the brightness observed along this path is generated by electrons 
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(a) 


Figure 3.5 a-c PERCENT OF THE TOTAL 
EMISSION COEFFICIENT J(f) 
(FOR GIVEN FREQUENCIES) AS 
A FUNCTION OF ELECTRON ENERGY 
FOR VARIOUS L VALUES. 


PERCENT OF THE TOTAL EMISSION COEFFICIENT PERCENT OP TOTAL EMISSION COEFFICIENT 



ELECTRON ENERGY (kev) 


(b) 



ELECTRON ENERGY (kev) 


(c) 
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Figure 3.6 a-c PERCENT OF THE TOTAL 
EMISSION COEFFICIENT J(f) 
(FOR GIVEN FREQUENCIES) AS 
A FUNCTION OF HARMONIC NUMBER 
FOR VARIOUS L VALUES. 


PERCENT OF THE TOTAL EMISSION COEFFICIENT PERCENT OF THE TOTAL EMISSION COEFFICENT 




HARMONIC NUMBER 


(c) 






with 100 kev < E < 500 kev and from harmonic numbers 3 <; s <20, For 
f = 300 kHz, 70% to 80% of the brightness is accounted for by electrons 
with 600 kev < E < 3 Mev and 50 < s < 600. 

This qualitative relationship between the electron energies, 
harmonic numbers, and frequency of observation was found to hold for 
each Kp range over the entire emission region. Consequently, if the 
gyro- synchrotron noise from the outer radiation belt can be observed, 
appropriate calculations should give an estimate of the density and 
energy spectrum (integrated over the emission region) of the electrons 
which radiate the power. For a study of this nature to be feasible, 
however, it must first be determined if the noise is sufficiently 
intense to be detected. Since this question is considered in Chapter V, 
further discussion of the results of this section will be postponed 
until that time. 

G. The Brightness Observed along Ray Paths off the Geomagnetic 
Equatorial Plane 

The calculations for the brightness observed along ray paths off 
the geomagnetic equatorial plane are more complicated than those of the 
preceding sections. The difficulties arise from the variation of the 
magnetic flux density (B) , observation angle (9), and the particle 
fluxes along segments of a ray path within a given L shell. (In the 
geomagnetic equatorial plane these quantities are constant within each 
L shell.) The analysis of this section can be simplified, however, by 
considering ray paths in the plane defined by the earth-observer line 
and the geomagnetic dipole axis (henceforth referred to as the 0-A 


112 



plane).' Limiting the rays to this plane reduces the problem from three 
dimensions to two dimensions, but still allows an accurate value for 
the average brightness of the magnetosphere (defined in Section III.G.3) 
to be calculated. 

1. Before calculating the average brightness it is 

necessary to obtain the angular extent of the emission region as a 

function of frequency and in the 0-A plane. This bound is determined 

in a manner similar to that used for deriving e (f) in Section 

max 

III.E.3. In this case, define a (f) as the angle which gives the ray 

max 

path in the O-A plane where = 0.5[l(a - 0° , f ) ] . (Note that 

the path given by a = 0° in this plane is the same as the one defined 
o 

by e = 0 in the geomagnetic equatorial plane.) Using this definition, 

the values of a (f) were calculated as a function of K and are given 
max P 

in Table 2 of Section III.G. 

The results for a (f) show that, as f increases, there is a sub* 
max' ’ ’ 

stantial decrease in the angular extent of the emission region in the 

OA plane, contrary to the results obtained for e (f). (Compare 

max 

and e (f) in Table 2.) This characteristic of a (f) is 
u max v 7 max' 7 max v ' 

explained by recalling the relationship established between the 

frequency of observation (f) and the range of harmonic numbers and 

electron energies which account for the radiation at f, To be more 

precise, Figure 1.10 illustrates that at low harmonic numbers, the 

power radiated by electrons with small pitch angles is appreciable 

(i.e., within a factor of one-half) compared to the power radiated by 

an electron of the same energy with o' = 90°. Furthermore, the magnetic 

flux density (B) increases (Equation 2.7) as \\ I increases from zero. 

1 m 1 
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Thus, tfce radio noise observed at low frequencies along ray paths in 
the O-A plane comes from harmonic numbers lower than those which 
account for the noise emanating from the geomagnetic equatorial plane. 
(The radiation observed at low frequencies along ray paths in the O-A 
plane results from very low harmonic numbers (1 £ s < 10) and hence is 
generated predominantly by the cyclotron emission process.) 

For these low harmonic numbers, Figures 1.10 and 1.7-1. 9 show that 
appreciable power can be radiated by electrons with 10° £ <y ^ 90° at 
angles 0 close to zero degrees. Thus, because low frequencies are 
generated at low harmonic numbers , appreciable power at these frequen- 
cies will be radiated along ray paths which pass nearly tangential to 
the field lines 5 £ L £ 7 in the O-A plane. In particular, the 
brightness observed at f a 30 kHz along ray paths passing tangential 

to the shells 6.5 < L < 7 (i.e., a (f=30 kHz) ~ 7°) is within 50% of 

max 

the brightness observed along the path a = e = 0° at this frequency. 

For the higher frequencies, however, the situation changes dras- 
tically. If we consider radiation at f = 300 kHz, the noise is 
emitted from high energy electrons, E> 600 kev , at harmonic numbers 
close to (i.e., mainly by a synchrotron emission process), where 

s pk * s < ^ e ^ : *- nec ^ Section I.D.l. As illustrated in Figure 1.10, the 

power radiated by a particle with a = 90°, at s = s (a = 90°), is much 

pk 

greater than that radiated by a particle of the same energy at the same 
harmonic number (and hence frequency), but with a < 70° . Furthermore, 
the beamwidth of the radiation at the higher frequencies is much less 
than that of the lower harmonics (see Figures 1.7 - 1,16). Therefore, 
contrary to the results for f = 30 kHz, appreciable power at f = 300 kHz 
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is obtained only along ray paths which make angles less than 30° with the 

normals to the field line L - 5 .0 in the 0-A plane. (This geometry 

results in electrons, with 70° < a < 90° , emitting radiation directly 

toward the observer along these paths.) For a spacecraft located at 
5 

2.04 x 10 km from earth in the geomagnetic equatorial plane this gives 

a a (f = 300 kHz) of 2,0°. The consistent decrease in a (f) as f 
max max 

increases is, then, a result of the emission process going from cyclo- 
tron at low frequencies to gyro- synchrotron at intermediate frequencies 
and, finally, to synchrotron at high frequencies. 

2. The calculations for CJ max (f) and for 1(f) observed along ray 
paths in the O-A plane showed that the line integral of Equation (3.2) 
could be limited to segments of the ray paths located in the portion 
of the magnetosphere bounded by \\ \ < A0 ° . Truncating the paths in 
this manner simplified the analysis of this section, while introducing 
at most a 10% decrease in the results obtained for I Qv (f) in III.G.3, 
an error which was felt to be insignificant. 

3. Using the results of preceding sections, it is now possible to 

consider the effects of the bulge in the plasmapause on the observed 
brightness. To this end, define the average brightness of the magne- 
tosphere as [Kraus, 1966] 

d V 


I (f ) = 

av v ' 


—fh- 

f ,K n ) JJ 


e> d0 


^ '“p 


(3.3) 



and 



(2a) (2e) 
2282 


(3.4) 


where 

I(a,e) = brightness observed along the ray path given by angles 
a and e . 

dQ = element of solid angle subtended by the source at the 
observer . 

0 (f ,K ) = solid angle subtended by the source at the observer, 
s p 

The quantity, 3282, converts degrees squared to steradians . 
The factors of two included with e and cr account for the symmetry of 
the source region (when the bulge is neglected) in both the 0-A plane 
and the geomagnetic equatorial plane. Using these equations, I av (f) 
was calculated as a function of K , and the results are given in 
Table 2. 

Since the bulge in the plasmapause is much more evident during 

geomagnetically disturbed periods (i.e., when the plasmapause, in the 

local time zone 2100 to 1800 hr, is compressed to L < 5), the 

discussion in this section shall be focused primarily on times when 

5 <; K £ 6 + . For this K range, the bulge will have its maximum effect 
P P 

on the observed brightness when it is located on the earth* s limb (as 
seen by the observer in the interplanetary medium). That is, at this 
location the bulge will shield the largest fraction of the emission 
region from the observer. 

To determine the fraction which cannot emit radiation toward the 
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Table 2 


SUMMARY OF RESULTS OBTAINED IN SECTIONS III.E.3, G.l, AND G.2 


Frequency 

(kHz) 

K p Range 

e (f) 

max 

(deg) 

a (f) 

max 

(deg) 

fi (f ,K ) 
s p 

(ster) 

I (f) 
a £ 

(watts/m -Hz-ster) 






„ -21 


lM 

VI 

M 

VI 

rH 

11.8 

7.0 

0.100 

5.5 X 10 


P 




-20 

30 

3 5 K 5 4 

11 .4 

6.9 

1.096 

6.5 X 10 


P + 




„ -19 


5 5 K 5 6 

10.8 

6 .7 

0.090 

5.3 X 10 


P 










, -21 


1 5 K 5 2 

11.6 

4.5 

0.064 

2.6 X 10 


P 




-20 

44 

3 5 K 5 4 

11.3 

4,4 

0.061 

3.1 X 10 


P + 




-19 


5 5 K 5 6 + 

10.7 

4.3 

0.056 

2.5 X 10 


P 










-21 


1 5 K 5 2 

11.4 

3.5 

0.049 

2 1,2 X 10 


P 




-20 

55 

3 5 K 5 4 

11.2 

3.5 

0.049 

2.3 X 10 


P + 




„ -19 


5 5 K 5 6 

10.7 

3.2 

0.042 

2.1 X 10 


P 










-21 


1 5 K 5 2 

11.3 

3.3 

0.045 

1.6 X 10 


P 




. -20 

67 

3 5 K 5 4 

11.0 

3.2 

0.043 

1.8 X 10 


P + 




, -19 


5 £ K ^6 

10.7 

3.1 

0.041 

1,3 X 10 


P 










, -22 


1~ 5 K 5 2 

11 .1 

2.9 

0.039 

9.4 X 10 


P 




-20 

83 

3 5 K 5 4 

10.9 

2.7 

0.036 

1.1 X 10 


P + 




_ -20 


5 5 K 5 6 

10.8 

2.5 

0.033 

7.7 X 10 


P 










„ -22 


l” 5 K 5 2 

11.0 

2.7 

0.036 

8,3 X 10 


P 




_ -20 

92 

3 5 K 5 4 

10.8 

2.5 

0.033 

1.0 X 10 


P . 




. -20 


5 5 K 5 6 

10.8 

2.4 

0.032 

6.6 X 10 


P 










„ „ , -22 


03 

VI 

W 

VI 

rH 

10.8 

2.5 ' 

0.033 

6.5 X 10 


P 




„ -21 

110 

3 5 K 5 4 

10.6 

2.3 

0.030 

8.3 X 10 


P . 




„ -20 


5 5 K 5 6 

10.7 

2.2 

0.029 

5.4 X 10 


P 

1 
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Table 2 


CONTINUED 


Frequency 

K Range 
P 

1 

8 

a 

Cl 

s 

I 

av 


l“ 5 K £ 2 

10.7 

2.:3 

0.030 

4.6 X 10" 22 

150 

P 

3 £ K 5 4 

10.5 

2.2 

0.028 

-21 

5 .5 X 10 


5 £ K < 6 

10.6 

2.1 

0.027 

3.8 X 10" 20 


p 

1 5 K 5 2 

10.5 

2.1 

0.028 

3.2 X 10" 22 

200 

U 

VI 

CO 

10.4 

2.0 

0.025 

-21 

4.0 X 10 


F + 
5 5 K 5 6 

10.7 

1 .8 

0.023 

3.1 X 10 -20 


P 

l” 5 K ^2 

10.5 

2.1 

0.027 

2.9 X 10" 22 

250 

P 

3 ^ K 5 4 

10.3 

1.9 

0.024 

3. 7 X 10 21 


P . 

5 5 K 5 6 

10.6 

1.8 

0.023 

2.6 X 10" 20 


P 

1~ < K ^2 

10.5 

2.1 

0.027 i 

-22 

2.6 X 10 

300 

p 

3 5 K 5 4 

10.3 

2.0 

0.025 

3.0 X 10 21 


P + 

5 5 K 5 6 

10.6 

i 

1.8 

0.023 

2.5 X 10 20 


P 




observer, we note that the bulge extends to L ^ 5.0 in the 1800 to 
2100 hr local time zone (see Section II. C). However, if we assume as 
a first approximation that it extends to L - 7.0 in this time zone, 
it is evident that the power radiated in one-half of the solid angle 
Q (f,K ) would be reduced to zero. Consequently, the average brightness 
would be decreased by 50%. 

To improve on this first estimate for the decrease in I (f). it 

av 

should be noted that at least 30% of the power radiated in the geo- 

+ 


magnetic equatorial plane (when 5 5 K ^ 6 ) comes from the region 

P 

5 £ L £ 7 (see Figures 3.4 a-c) . Thus the decrease in I (f) will, in 

a v 

reality, be no more than 35% (i.e., 1.9 dB). Since this decrease in 
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the average brightness is well within the approximations used through- 
out the analysis presented thus far, it is concluded that the bulge 
in the plasmapause can be neglected in this study. 


H. The Average Calculated Spectrums of the Radio Noise Generated in 

the Earth* s Outer Radiation Belt 

The values for I ( f ) given in the last column of Table 2 have 
av 

been plotted in Figures 3,7 a-c to illustrate the spectrums calculated 
for the gyro- synchrotron radio noise. Two important aspects of these 
figures are discussed below. 

1. Each spectrum can be approximated by two straight lines, one 
for 30 kHz £ f £ 92 kHz and one for 110 kHz £ f £ 300 kHz, which differ 
by a factor of approximately two in their slopes. In Section III.G.l, 
we showed that, for a frequency (f) , the noise is emitted by electrons 
with ^ o' £ 90°, where decreases monotonically as f increases. 

(Recall that ^ 10° for f : 30 kHz and increases to > 70° for 
f > 110 kHz.) Since the number of particles of energy E radiating at 
a given frequency along the direction of observation is given approxi- 
mately by (Section I.F) 


N(E) = 2tt 


i 


N(E,or) sin a dev 


(3.5) 


N(E) also decreases as (f) increases. The magnitude of |d[log N(E)]/ 
d[log f]|, however, is less for f £ 110 kHz than for 30 kHz £ f £ 110 
kHz (i.e., the emission process changes rapidly from cyclotron 
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AVERAGE CALCULATED BRIGHTNESS (wotls/m 2 -Hz-sfer> 



FREQUENCY (kHz) 


00 (b) (c) 

Figure 3.7 a-c THEORETICAL SPECTRUMS FOR THE GYRO- SYNCHROTRON RADIO NOISE GENERATED IN THE 
EARTH'S OUTER RADIATION BELT. In this figure (a) is the exponent in 1(f) « f 





(f < 40 kHz) to gyro- synchrotron (44 kHz < f < 92 kHz) but becomes, 
and remains, synchrotron for f > 110 kHz.) Hence, for 30 kHz £ f £ 92 
kHz, |d[log I av (f)]/d[log f]| is necessarily greater than for f > 150 
kHz. (The value 0.9 obtained for | d[log I Qv (f ) ]/d[log f]| when 110 
kHz < f < 300 kHz agrees very well with the value calculated by Vesecky 
[1969], ~1.0, which he obtained for the synchrotron noise spectrum, 

5 MHz < f < 20 MHz, generated by electrons trapped above the earth 1 s 
auroral zones.) 

2. The spectrums in these figures are also characterized by only 

a slight increase in ]d[log I qv (f ) ]/d[log f]| as increases from 

1 ^ Kp £ 2 to 5 < Kp < 6 + # even though the average brightness for 

30 kHz < f < 300 kHz increases by almost two orders of magnitude. This 

feature of the data is to be expected, since the electron fluxes in the 

outer radiation belt are enhanced nearly uniformly over all energies 

during sustained geomagnetic activity or during the early recovery 

phase of large geomagnetic storms (Section II. B). (The slight increase 

in | d[log I av (f )]/d[log fjj is accounted for by more low energy 

electrons radiating appreciable power at low frequencies when 

increases. That is, regions of higher magnetic flux density, and hence 

f^, are exposed as the geomagnetic activity increases. Consequently, 

noise at low frequencies (f^) i s augmented by radiation from electrons 

for which f . « These particles do not radiate appreciable 

Jo pK e 

energy at for lower K p values, since then f^ = sf^ where s > s ; 
see Figure 1 .6) . 

There exists one possibility (neglected to this point) which could 
affect the constancy of |d[log I gv (f )]/d[log f]|: deformation of the 
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dipole field geometry due to the compression of the magnetosphere by 

the solar wind. During severe disturbances, the dipole field lines 

are compressed at the subsolar point to L «6 [Roederer, 1970]. This 

study, however, has considered the portion of the magnetosphere within 

the L « 7,0 shell and moderate geomagnetic disturbances resulting in 

5 £ K £ 6 + , It is therefore expected that the dipole field ’geometry', 

assumed for L < 7, is a valid representation of the earth’s magnetic 

field. Thus, the slight increase in |d[log I gv (f )]/d[log f]| for 

increasing K is felt to be a valid characteristic of the gyro- 
P 

synchrotron radio noise spectrums. 
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CHAPTER IV 


EXPERIMENTAL DATA 


A . The Spacecraft 

The IMP-6 (Explorer 43) was launched March 13, 1971 into a highly 
elliptical orbit defined by an eccentricity of 0.94, an initial perigee 
of 354 km, and an initial apogee of 206,000 km. (32 earth radii). The 
duration of this orbit was 4.18 days and it was given an inclination 
of 28.7 degrees. 

The spacecraft is a sixteen-sided drum which is spin stabilized 
about the drum axis at a rate of 5.4 revolutions per minute. The spin 
axis is perpendicular to the ecliptic plane and points toward the 
southern ecliptic pole. 

B. The Radio Astronomy Experiment 

The GSFC radio astronomy experiment aboard the IMP-6 is composed 
of two total power, stepped frequency radiometers. Each receiver has 
thirty- two discrete frequency channels which are sampled over a 5.11 
second interval. They are not, however, entirely similar. One covers 
the frequency range 30 kHz to 9900 kHz, has a bandwidth of 10 kHz, and a 
time constant of 6 milliseconds. The second has a frequency range 
from 30 kHz to 4900 kHz, a bandwidth of 3 kHz, and a time constant of 
40 milliseconds. 

Both receivers are connected to the same 300 foot dipole antenna, 
which is situated in the spacecraft spin plane. Due to this location 
and the satellite spin, any source of radio emission not located at 
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one of the ecliptic poles will have the maximum or minimum of the dipole 
antenna pattern swept past it, resulting in a modulation of the received 
signal. By analyzing the depth and phase of this modulation, the 
apparent size and location of the source can be determined [Brown, 1972], 

To measure the temporal stability of the radio receivers, an 
automatic internal calibrating system was included in the radio astron- 
omy package. The data returned from this system are used to compensate 
for time variations in the radio noise measurements due to short and 
long term drifts in receiver gains. 

C. Initial Data Reduction 

Dr. Larry Brown [1972] has analyzed the raw data from the tenth 
through the fortieth orbits of the IMP-6 to determine the steady 
underlying radiation observed by the spacecraft. His data reduction 
technique was specifically intended to reduce the contamination of the 
background radio spectrum by sporadic emissions from the sun, earth 
and the spacecraft itself. His results showed that three distinct 
forms of background radiation could be identified. 

The first form of emission was that generated by galactic and 
extragalactic sources. This component, termed the radio background 
spectrum, was studied in detail by Brown [1972] for frequencies between 
200 kHz and 1700 kHz . 

The second form of radiation, found to emanate from the earth, was 
sporadic in nature and showed a peak in its brightness spectrum at 
about 250 kHz. This spectrum showed a decrease by two orders of 
magnitude within tens of kilohertz on either side of the maximum. 
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The third was a continuous, highly variable low frequency compo- 
nent, observed between 30 kHz and 110 kHz. Over this frequency range, 
the average spectrum for this radiation satisfied an expression of 
the form 

I (f) = 3.3 x 10" 8 ±~ 2 S (4.1) 

m 

where 

-2 -1 -1 

I = measured source brightness in watts m Hz ster 
m 

f = frequency in Hz, 

This result for X was obtained by assuming the source filled the IMP-6 
m 

dipole antenna pattern [Brown, personal communication]. Furthermore, 
the direction finding capabilities of the radio astronomy experiment 
showed this radio noise also to come from the earth. Spectrums 
representative of these three forms of background radiation are shown 
in Figure 4.1. 

In this chapter we shall derive a number of important features 
of this third type of emission. These characteristics were determined 
by analyzing hourly plots of the background radiation observed by the 
IMP-6 during May through October 1971 (for which the author expresses 
his gratitude to Dr. Brown who supplied the data in this reduced form) . 
However, the choice of plots used in this paper was restricted to times 
when the spacecraft was at least 194,000 km from the earth. This 
restriction was to provide an observer geometry similar to the one 
used in the theoretical analysis of the previous chapters. The amount 
of IMP-6 data left out is not large, as the spacecraft spends a large 
fraction of its time in each orbit at these distances from the earth. 
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Figure 4.1 BACKGROUND RADIATION OBSERVED BY THE IMP-6 DURING APRIL- 
AUGUST 1971. (THESE DATA WERE CALCULATED FROM IMP-6 
TEMPERATURE MEASUREMENTS ASSUMING THE GALACTIC SOURCES 
AND THE SOURCES LOCATED AT THE EARTH FILLED THE SPACECRAFT 
DIPOLE ANTENNA PATTERN.) [This figure was supplied by 
R. G. Stone, personal communication, 1972.] 
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(This can be proven by using Kepler’s law of planetary motion [Bate et 
al . , 1971] and the orbit described in Section IV. A). 

D. Correlation between the Radio Noise Generated at the Earth and 
Geomagnetic Activity 

The most salient feature evident in the data gathered by the 
IMP-6 was a two order of magnitude variation in the brightness for the 
radiation between 30 kHz and 110 kHz. To determine the cause for this 
variation in noise intensity, the IMP-6 data were correlated with geo- 
physical and solar phenomena. It was found that an increase in the 
intensity of the radiation generally followed the onset of geomagnetic 
activity, and as the disturbance subsided the radiation returned to its 
unperturbed value. To put this observation into more precise terms, 
average spectrums for the noise as a function of were calculated by 
subdividing the data into three groups. The first and second contained 
data for periods of sustained geomagnetic activity given by 1~ £ £ 2 

and 3 £ Kp £ 4 respectively. The third group comprise, data taken 
during the early recovery phase of the 16 - 18 May and 24 - 28 September 
1971 geomagnetic storms. (The storm characteristics for the first case 
are shown in Figure 4.2). During this phase of the storms, the index 

was bounded by 5 £ K £ 6 + , (For the May storm, K reached values as 

P P 

high as 7 for a very short period of time. The radio noise data for 

this high value of K^are, therefore, rather sparse and will be used 

only for argumentative purposes in the next chapter.) The resulting 

average brightness spectrums (as functions of K ) are given in 

P 

Figures 4.3 a-c. 
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MAGNETOGRAMS of geomagnetic storms 


Figure 4.2 MORPHOLOGY OF THE GEOMAGNETIC STORM OCCURRING ON 
16-18 MAY 1971 [From Solar Geoph 3 ^sical Data]. 
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BRIGHTNESS (wotts/m -Hz-ster) 



Figure 4.3 a-c BRIGHTNESS SPECTRUMS FOR THE SOURCE LOCATED AT THE EARTH:. (THESE DATA WERE 
CALCULATED FROM IMP-6 TEMPERATURE MEASUREMENTS ASSUMING THE SOURCE FILLED THE 
SPACECRAFT DIPOLE ANTENNA PATTERN.) In this figure (a) is the exponent 1(f) « f~ 





It should be noted that in Figure 4.3a, the experimental data were 
extrapolated for f > 55 kHz. This was necessary as the radio noise 
power for frequencies above 55 kHz was of such small amplitude that the 
receivers w r ere operating near the lower limit of their dynamic range. 

At these power levels, the linearity of the receivers is in question 
[Brown, personal communication]; therefore, the data for f -» 55 kHz were 
not used* Nevertheless, there were sufficient data for f =30 kHz to 
55 kHz to obtain the curve shown in Figure 4.3a. 

These figures illustrate that the spectrums for each range 
follow a law of the form 

I (f) = I f" a (4.2) 

av c 

where (a) increases as the K value increases. If the entire set of 

P 

orbital data covering April through August is averaged, the resulting 

spectrum is indeed given by Equation (4.1). 

The error bars in these figures represent one standard deviation 

in the data about the average points plotted. Thus, even though the 

spectrums correlate with changes in K values, there still exists a 

P 

large variation within each K range. These variations can be attri~ 

P 

buted to fluctuations of the trapped particle populations due to changes 
in the solar wind density and velocity as well as other geomagnetic 
related phenomena [Vette, 1971], 

For the lower K range, the fluctuations in brightness due to these 
P 

mechanisms should be a minimum. The error bars in Figure 4.3a, show 
however, that changes by a factor of 2.5 in the observed brightness are 
common. Even though this factor is less than those obtained for the 
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other two K ranges, it is still appreciable. Therefore, the spectrum 
P 

for 1 £ K ^ 2 will be reconsidered in Chapter V. 

P 

E. Determination of the Solid Angle Subtended by the Source Observed 
by the IMP-6 

In the next chapter, we shall compare the average theoretical 

brightness of the magnetosphere I (f), as calculated in Chapter III, 

a v 

to the average brightness of the earth (between 30 kHz and 110 kHz) 
measured on IMP-6 , It must therefore be determined if the source 
observed by the spacecraft fills the satellite dipole antenna* 
pattern, as assumed by Dr. Brown. This is necessary since the actual 
average brightness of the source I (source) is given by [Kraus,. 

a Cl 

1966] 

CJ (f ,K ) < n (f) (4.3) 

S p 

n s (f ,K p ) > n A (f> (4.4> 


I (source) = 
act 


fi(f ) 


n s (f 'V 


I (source) 
m 


or 


I (source) = I (source) 
act m 


where 

I (source) = actual measured brightness of the source obtained 
act 

. from IMP-6 experiment. 

I^(source)= average source brightness calculated by assuming 
the source fills the spacecraft antenna beam. 
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Q (f k ) = the solid angle subtended by the source at the observa- 
s ' p 

tion point for a frequency (f) and given K value. 

ir 

0 A (f) = the solid angle of the spacecraft antenna pattern 

This question was resolved by plotting the IMP-6 data for f = 30 kHz, 

44 kHz, and 55 kHz (when S 3) as a function of distance from the 

earth to the satellite. (A lower limit of ten earth radii was used in 

this analysis to avoid problems related to the spacecraft being within 

the earth's magnetoplasma . ) The results are shown in Figures 4.4 a-c. 

If the source subtends a solid angle less than that of the antenna 

pattern, the brightness should decrease by about an order of magnitude 

as the distance increases from ten to thirty-two earth radii (i.e., the 
_ o 

R decrease in the radio flux density should be observable). The 

results for f = 30 kHz (Figure 4.4a) are inconclusive due to the 

scattered nature of the IMP-6 data. (The fluctuation in these data will 

be discussed in the next chapter.) For f - 44 kHz and 55 kHz, however, 

the corresponding figures do show a decrease in the measured brightness 

by approximately one order of magnitude as the spacecraft distance 

from the earth increases by a factor of 3.2. Since the same results 

were also obtained at the higher frequencies, it must be concluded that 

the source actually does not fill the satellite antenna pattern when 

4 

the spacecraft is located more than 6,4 x 10 km from earth. Thus, if 
a valid comparison between the IMP-6 data and the theoretical calcula- 
tions is to be made, either set of data must be corrected. 

Since the calculated data give the actual brightness of the earth’s 
magnetosphere, while the IMP-6 data were calculated on an assumption, 
this latter set of data shall be corrected. The correction factor, 
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Figure 4.4 a-c SCATTER PLOT OF 
THE BRIGHTNESS OF THE 
EARTH (CALCULATED FROM 
IMP-6 DATA UNDER THE 
ASSUMPTION STATED IN 
THE CAPTION OF FIGURE 
4.3) AS A FUNCTION OF 
THE SPACECRAFT GEOCEN- 
TRIC DISTANCE. 






0 (fVQ (f K ) is obtained by assuming that the IMP-6 data are in fact a 
A s ' p 

measure of the gyro- synchrotron radiation emanating from the earth's 
outer radiation belt (i.e., 0 (f ,K ) will be taken from Table 2). Then, 

S p 

if the corrected data agree with the theoretical calculations, the 
correction is justified and the recalculated spacecraft data will also 
give the actual brightness of the earth’s magnetosphere. To this end, 
the theoretical results of Figures 3.7 a-c will be compared to 
Figures 4.5 a-c, which show the corrected IMP-6 data. 
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FREQUENCY (kHz) 

Figure 4.5 a-c CORRECTED IMP-6 BRIGHTNESS SPECTRUMS (COMPARE TO FIGURES 4.3 a-c) WHEN THE SOURCE 
SOLID ANGLE, f2 s (f,Kp) IS ACCOUNTED FOR. In this figure (a) is the exponent in i(f)ocf“ 





CHAPTER V 


COMPARISONS OF THEORETICAL AND EXPERIMENTAL RESULTS, 

FINAL CONCLUSIONS AND APPLICATIONS 

A Comparisons of Results 

Equation (1.15) indicates that the brightness observed along ray 
paths passing through the earth’s magnetosphere is a result of the 
radiation emanating from within this region and from galactic and extra- 
galactic sources. The IMP-6 data (Figure 4.1) show, however, that the 
emission from the latter sources for 30 kHz £ f £ 110 kHz is orders of 
magnitude below that from the earth's environment; consequently, the 
galactic and extragalact ic radio noise can be neglected in the discus- 
sion to follow. 

There exists one problem pertaining to observer locations which 
must be considered before a comparison can be made between the theoret- 
ical results and the experimental data. In Chapter III the calculations 
were based on an observer situated in the geomagnetic equatorial plane. 
The IMP-6 orbit, however, is inclined nearly 30° to the earth's geo- 
graphic equatorial plane. Consequently, the spacecraft’s inclination 
angle (0 inc ) the geomagnetic equatorial plane varies diurnally from 
«s20° to « 40°. 

Thorne [l963] conducted a detailed study on the variation expected 

in the synchrotron radiation from stars as a function of an observer’s 

inclination angle to their magnetic equatorial planes. His results 

show that, for a magnetic field geometry identical to the one considered 

in the paper, the noise intensity is maximum in the equatorial plane 

(0 . = 0°)and decreases monotonically to a minimum value for 0 = 90° 

me inc 
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He found, furthermore, that ] I av © ±nc = 90 °)/ I av ( e lnc = °°>| is a func_ 

q 

tion of the exponent (q) in the sin a pitch angle distribution he 

assumed for the star's stably trapped electrons. In particular, when 

0 increased from zero to forty degrees, the noise intensity decreased 
inc 

by a factor of two for q = 6.0 but only by a factor of 1.1 for q = 2.5. 

2 

In this study, a sin a pitch angle distribution was found appro- 
priate for the particles trapped in the earth's magnetosphere (when 3 < 
L < 7) . Furthermore, I dl (f )/d0 . I is expected to be greater for 
synchrotron radiation than for cyclotron radiation. (For cyclotron 
emission the particles radiate nearly isotropically, while for synchro- 
tron emission they radiate predominantly along the normals to the 
magnetic field lines, see Sections I.D and III.G.) Therefore, Thorne's 
factor of 1.1 can be considered an upper bound for the discrepancy 
expected between the theoretical results and the experimental data due 
to the different observer locations. Since this factor is insignifi- 
cant, a direct comparison of Figures 3.7 a-c and 4.5 a-c (which are 
redrawn in Figure 5.1 for convenience) can be made. A summary and 
discussion of the results obtained from this analysis are given below, 

1 . For 30 kHz £ f <, 110 kHz the slopes for the theoretical noise 

spectrums in Figures 3.7 a-c were averaged for the three K ranges, and 

P 

the value ]d[log I (f ) ]/d[log f]|^ e °was then compared to a similar 
average obtained for the corrected IMP-6 spectrums in Figures 4.5 a-c. 
The result obtained was 


K log * av ( f >] /d [i°e f ]| av 


(1.05) |d[log I av ( f )]/d[log f] 


theo 

av 
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BRIGHTNESS lwott$/m - Hz ~ ster ) 



Figure 5.1 a-c AVERAGE THEORETICAL AND EXPERIMENTAL SPECTRUMS FOR THE RADIO NOISE EMANATING FROM 
THE EARTH’S OUTER RADIATION BELT. In this figure (a) is the exponent in Iff) <x f a 





Furthermore 4 when the slopes for each K range were compared separately, 

P 

the largest discrepancy occurred for 5 <, £ 6 + where 

| d[log I av (f)]/d[log f]^ IP ~ 6 = (1.09) |d[log I av (f)]/d[log f]| a J e ° 

2. Both sets of data show a slight, but discernable, increase in 

|d[log I (f)]/d[log f]| as K increases, 
a v p 

3. The increase in the magnitude of I^(f) in the IMP-6 data 

(which was correlated with increases in K ) also exists in the theo- 

P 

retical calculations. In addition, I dl (f)/dK I ® is within 5% of 

1 av p 1 

I dl (f)/dK I " theo for each K range. Thus, the average enhancement of 
1 av p 1 p 

the stably trapped electrons due to geomagnetic disturbances satisfac- 
torily accounts for the increase in the radio noise observed by the 
IMP-6 during such activity. 

4. In the worst case, I (f ,K ) = (15.0) I (f,K ) 

’ av * p IMF-6 av v ' p theo. 

This constant discrepancy, however, is not felt to be critical for the 
following reason s. 

i) Throughout this paper, approximations, worst case condi- 
tions, and assumptions have been used to simplify the analysis which led 

to an underestimation of I (f ) . Individually, each of these simplifi- 

a v 

cations results is no more than 10% decreases in I (f) . Their combined 

av 

effect, however, results in a total reduction of I (f) by approximately 

av 

a factor of two (due to the number of simplifications made). 

ii) In Chapter III, calculations for the ordinary mode were 
not presented since the power radiated into it is generally less than 
or equal to that radiated into the extraordinary mode. However, 
calculations for free-space showed that the power radiated into the 
ordinary mode is comparable to that radiated into the extraordinary mode 
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when 1 9 | , > [ ot - 10°| . Since the emission region is nearly free-space 

for 30 kHz £ f £ 110 kHz, the power radiated into the ordinary mode will 

increase the values calculated for I (f ,K ) by about 1.5. 

av p 

iii) in this study, superthermal electrons in the pseudo- 
trapping regions [Vette, 1971] have been entirely neglected due to the 
lack of sufficient data on the particle distributions throughout this 
segment of the magnetosphere. However, measurements have been made for 
the fluxes in the portion of the pseudo-trapping region given by 8 £ L 
£ 10 in the morning magnetosphere [Hess, 1968; Williams and Palmer, 
1964], The data show that at high geomagnetic latitudes (and hence 
regions of high magnetic flux density) the integrated electron fluxes 
in this region (for E > 40 kev and E s 280 kev) are two orders of 
magnitude greater than the fluxes in the corresponding segment of the 
nighttime magnetosphere. Thus, fluxes in the former region are 
comparable to those which exist at much lower L shells in the latter 
(i.e., 6 < L < 7) . The pseudo -trapping region in the morning magneto- 
sphere should therefore affect the results of this study in two ways. 
First, the radiation from electrons in this region will increase the 

values calculated for I (f,K ). Second, the symmetry of the emission 

av p 

region will be destroyed (see Section III.G.2) since a larger segment 
of the morning magnetosphere will contribute detectable radiation, and 
the electron fluxes in 8 < L < 10 will make its overall (spatially 
integrated) brightness considerably greater than that of the nighttime 
magnetosphere. Due to lack of data, however, quantitative values for 
these effects can not be calculated. 

From the above discussion, it is evident that the calculations for 
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I (f ,K ) can be below the actual brightness of the magnetosphere by 
av p 

over a factor of three. Hence, the noise detected by the IMP-6 and the 
intensity of the radiation emanating from the earth 1 s magnetosphere 
differ not by a factor of 15, but more likely by a factor of less than 
five, i.e., a discrepancy of at most 7 dB in the noise power, 

B Summary and Conclusions 

In this study the problem of detecting gyro- synchrotron radio 
noise, from superthermal electrons trapped in the earth’s magneto- 
sphere, has been investigated for the frequency range 30 kHz to 300 kHz. 
The energy radiated by these electrons and propagated to an observer 
in the interplanetary medium was calculated using the equation of 
radiative transfer. To solve this equation, models for the thermal 
and superthermal electron densities as a function of geomagnetic 
activity were derived from in situ and indirect measurements of 
these particle populations. 

The analysis shows that appreciable radio noise in the LF range 
is emitted by electrons in the earth’s outer radiation belt. The inner 
belt, however, does not contribute power at these frequencies since it 
is located in the plasmasphere , a region of high thermal electron 
density which inhibits radiation at these low frequencies. The inten- 
sity of the radio noise from the outer belt was found to be positively 
correlated with K^. This correlation is due to the enhancement of 
stably trapped electrons during periods of sustained geomagnetic 
activity. For a frequency f in Hz (30 kHz <. f £ 150 kHz), the average 
brightness I (f) is given by (see Figure 5,1) 
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(5.1) 


-a -2 -1 

I (f) = M f watts m Hz ster 

av 

where both M and (a) vary with as follows: 

l’ <: K £ 2 
P 

3 £ K £4 
P 

5 <; K 5 6 + 
p 

The theoretical results were compared to data obtained from the 
GSFC experiment aboard the IMP-6 (as shown in Figure 5.1 and discussed 
in the preceding section) . The models were found to accurately predict 
the spectral shape, temporal variations and magnitude of the radio 
noise observed by the radio astronomy experiment. This author, there- 
fore, concludes that the electrons, which are injected and trapped by 
natural process in the earth’s outer radiation belt, generate detectable 
gyro-synchrotron radio noise in the LF range. The intensity of this 
noise is a function of geomagnetic activity, increasing as increases. 
Furthermore, this noise is observed only in the LF range since for 
f £ 2 MHz its intensity is below the cosmic noise level (compare Figure 
4.1 with Figure 5,1). 

0. Applications 

The gyro- synchrotron radiation at 110 kHz ^ f ^ 300 kHz was shown 
to come from high energy electrons, while the radiation at low frequen- 
cies, 50 kHz £ f £ 53 kHz, results from low energy electrons (Section 
III F) . During the main phase of a geomagnetic storm, the low energy 
electron fluxes, 40 kev £ E ^ 200 kev, are enhanced over quiet time 


M = 1 .51 x 

10 

a = 1.66 

M = 3 .16 x 

i — 1 

1 

o 

1— { 

a = 1.71 

M = 7.41 x 

io 11 

a = 1.82 
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values, while the high energy electron fluxes, E> 800 kev, are deple- 
ted (Section II. B) . Consequently, during this phase of a geomagnetic 
disturbance the gyro- synchrotron noise should have a spectrum whose 
slope is greater than those obtained in Figures 4.5 a-c. This is the 
case, as illustrated in Figure 5.2, which shows IMP-6 data taken during 
the main phase of the 16 - 18 May 1971 geomagnetic storm. 



Figure 5.2 CORRECT IMP-6 BRIGHTNESS SPECTRUM FOR THE NOISE EMANATING 

FROM THE EARTH'S OUTER RADIATION BELT DURING THE MAIN PHASE 

OF THE 16-18 MAY 1971 GEOMAGNETIC STORM, In this figure (a) 

- a 

is the exponent in 1(f) a f 
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Based on the discussion of the last paragraph, which summarizes an 
important result of this study, it is evident that the IMP-6 noise data 
can be used to obtain the average (spatially integrated) flux spectrum 
of the energetic electrons trapped in the earth 1 s outer radiation belt. 

A study of this nature would be conducted by first averaging noise data 
over long periods of time, then reversing the analysis presented in this 
paper to obtain electron flux spectrums similar to those given in 
Figures 2.1 and 2.2b. Furthermore, this type of study would allow the 
average temporal variations in the electron fluxes to be determined. 

For example, temporal variations in the radio noise data are evident 
in Figures 4.3a and 4.4a. During geomagnet ically quiet days, these 
figures show large fluctuations in the low frequency radio noise, 
f = 30 kHz, compared to the variations in the high frequency noise- This 
difference between the magnitude of the fluctuations is expected, as 
satellite measurements of the electron fluxes show consistent, large 
temporal variations in the low energy electron fluxes in contrast to 
stable , non-varying high energy fluxes [Freeman, 1964]. 

A second application of the results of this study is the use of 

the gyro- synchrotron radio noise to monitor the interplanetary solar 

wind density in the vicinity of the earth. This type of study is 

possible, since noise generated by this process should extend to 

frequencies below 15 kHz, and the interplanetary electron density at 

-3 

1.0 AU is an average of w 5 electrons- cm (i.e, , f p* 20.1 kHz) 

P 

[Brandt, 1970] . Thus, the medium will inhibit the propagation of the 
low frequency gyro- synchrotron radiation; the cutoff frequency of the 
noise will then determine f , from which the electron density can be 
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calculated, (In this calculation, the observer must be located near the 
earth in order to detect the radio noise. Therefore, the effects of the 
interplanetary plasma on the receiving antenna must be considered. 
Nevertheless, appropriate calculations should give an acceptable esti- 
mate of the electron density.) 

The above line of reasoning is supported empirically by the radio 

noise data, given in Figure 5.2, which show an unexpected roll-over in 

the spectrum at f « 47 kHz. For the day in which these data were 

taken (17 May 1971), the interplanetary electron density showed the 

following temporal variation (as measured by experiments aboard the 

-3 

IMF^6) . Early on the 17th, the density was 10 electrons-cm , which 

-3 

increased to 25 electrons-cm near noon, and subsequently returned to 

_ 3 

m 10 electrons-cm later in the day. Furthermore, for the period of 

time during which the data of Figure 5.2 were taken (0538-0638 UT) , the 

_ 3 

electron density was ^ 18 electrons-cm [Dr. S. J, Same, personal 
communication, 1973], resulting in f^ & 38 kHz. Consequently, as 

^ fp, the decrease in noise intensity for f < 47 kHz was probably 
due to propagation limitations (and, in this case, effects on the 
spacecraft antenna) created by the interplanetary plasma. (This appli- 
cation, and the one given below, can be implemented independent of the 
mechanism generating the noise detected by the IMP-6 experiments.) 

A final application involves the ray paths calculated in Section 
II. H. These paths, for 30 kHz < f ^ 300 kHz, do not penetrate the 
plasmasphere . Thus, through an occultation experiment, two satellites 
could be used to determine the location of the plasmapause . In such an 
experiment, one spacecraft transmits signals at f ^ 30 kHz to the IMP-6 
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spacecraft . Then by choosing an orbit for the former spacecraft which 
allows it to be intermittently occulted from the IMP-6 by the earths 
magnetosphere, the positions of the spacecrafts at occultation and ray 
path geometry can be used to determine the location of the plasmapause. 

The details of these applications have not been presented, as they 
would constitute another entire study. The point is, the gyro- 
synchrotron emission mechanism does satisfactorily account for the 
radio noise detected by the IMP-6 in the frequency range kHz <• f < 110 
kHz. Furthermore, observations of this radiation and the results 
obtained in this paper are useful in the study of the electrons trapped 
in the earth’s outer radiation belt, as well as in the study of other 
geophysical phenomena. 
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